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APPENDIX A 
DATA REDUCTION PRINCIPLES 


1. PROBLEM FORMULATION 

The simultaneous reduction method with constraints is achieved by finding 
values for all variables which minimize a functional comprised of the weighted 
sum of squares and differences between measured data and the observation model, 
and differences between the solution vector and a vector of apriori information 
about all of the parameters being estimated. In general if the vector Q is 
comprised of all of the parameters being estimated and the observational equation 
is given by: 

Map Gap(Q) 


where May is the measured quantity at the ju ground station of the zth aircraft 


location. 

The function G(@) relates the parameters being estimated to the quantities 
being measured and, therefore, forms the observation equations for the system of 
data reduction. For the LRPNS the quantities being measured are range changes at 
each ground station location. The data from each ground station is then a sequence 
of time differences which represent the changes in propagation time for each read 
command received from the aircraft. By suitable scaling, these time differences 
may be considered as representing changes in range between the first received read 
command and any succeeding received read commands. 

For purposes of discussion, then, the data available to be reduced represents 
the range change between the first aircraft location seen and all successive air- 
craft locations seen. The ground receivers each contajn a local oscillator which 
is used as a standard to measure elapsed time and hence range change. Since the 
ground station oscillator is not at exactly the same frequency as the aircraft 
reference oscillator and may be drifting at a fixed rate in frequency, the measured 
data will contain the effects of these discrepencies. The observation equation 
then contains terms which reflect not only the frequency offset and drift, but also 


an arbitrary initial range (or phase) term which corresponds to the range from the 
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ground station to the first aircraft location seen. 


The initial range parameter for the jm ground station is designated Ap 7 
The effect of the frequency differences between the aircraft oscillator and the jm ] 
ground station oscillator is contained in a term designated ays The drift rate 
between the aircraft and ground station oscillator is designated ao, for the j* i} 
station. ae 
The observation equation Gey is then given by ~~ 
Geet, = gay Gees Ae eue ij 
ge sak “Og 1j°k 2-5 Hl 
where 
- - 
if c 2\ 1/2 t 
m=1 rf 
a, is the initial range parameter ~ 
ay; is the frequency offset parameter | 
a. is the frequency drift rate parameter Tt 
ty is the time between the first and the zth aircraft interrogation command, 
It is worth note at this point to indicate that in fact the effects of the o 
oscillator offset and drift will be dependent on the time difference between the 5} 


first aircraft interrogation time seen by a particular ground station and the suc- 
ceeding interrogations. Thus, a particular ground station might first see the 
aircraft on the second, third or any later interrogation time rather than the first. 
It would be possible to include this exact situation in the data reduction only at 
the expense of a considerable rearrangement of the structure of the matrices involved. =; 
In so far as the positions of the ground stations such a reorganization is unnecessary 
since an arbitrary shift in the zero time point only changes the values of Ayr By or j] 
and ay recovered by the data reduction procedure. Therefore, to preserve the matrix 
structure (without detrimental effect on accuracy) the values of Aner ay, and ao. will 
be estimated as though all ground stations saw the first aircraft interrogation 7 


(whether or not any of them actually did). The only implication of this procedure is l 


that if the actual values of a; and ay. are desired (for test purposes perhaps) then f] 
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the results of the data reduction procedure will require some additional interpretation 
beyond a direct readout of the estimated values. 

The apriori information about the estimated parameters is contained in a vector 
@, and a covariance matrix Ww, which contains the estimates of the error in the apriori 


vector Q,° The functional to be minimized by a choice of @ is then given by 

F = [4 - o(@)|7w,[w - ca] + fe, - a]? |e, - o] , 
where Wy is the covariance matrix of measurement error and "5 is the covariance 
matrix of the apriori errors in Q, The minimization of the functional F is 
accomplished by partially differentiating F with respect to each of the parameters 
being estimated then equating the resultant set of (nonlinear) equations to the zero 
vector and solving for the value of Q which satisfies these equations. Then 


a =a"s OLA [ - c(a)] -W, [@, hs al 


where A(Q) is the matrix of partial derivatives of the function G(@) with respect to 


the vector @ evaluated at the value of Q. At the solution point the vector @ may be 


represented as 
Q=Q, + 89 


where Q,, represents the true values of the parameters being estimated. The function 


G(@) is then 
G(Q) = G(Q_) + 39,50 = G(@p) + A(Qp) 60. 

Substituting this value for G(Q) into the above relationship gives 
0 =-A7(Qp) Wy [M - G(@p)- A(Q)a] - Wy (@, - Op - 60) . 


The vector M - G(Qn) is, by definition of the observational equation, the vector of 


measurement errors 6. The vector @, = Qp will be defined to be 89 for convenience 


(it is the vector of errors in the apriori vector Q,)+ Making these substitutions 


and rearranging the equation gives: 


[a7 (o,) W, A(Q,) + we =a (Qp) W,6R + W560. 


As3 


The solution of this equation for 6Q would give the errors in the estimated para- 


meters (i.e., 56Q) as a function of the measurement error, 6R, and the apriori 


information errors, &Q. This equation also provides a basis for solving for the value 


of the vector @ which satisfies the set of nonlinear equations derived above. If 


6R is interpreted as the current difference vector between the measurements and the 


current value of the observation equation, and 6Q is the difference vector between the. 


apriori estimate of @ and the current value of @ a recursive method of finding Q 
can be formated. Suppose a, is the present value of @ then 


OR, =M- G(@,) A 


Then 
T A A be T LN be 

[4 (G5) W, AC Ww ]0, = 47(0,) WSR, + W780, , 
defines a value for 5Q, which may be used to determine a better estimate of Q; or 
a Rd 

The overall procedure then to start with ar /.it.al estimate of Gg 722 
(the apriori information), the two covariance nai vices Hi" and Ws and then 
iterate the value of the vector Q; until the correction 6Q, is the zero vector. 


The matrix of partial derivatives of G(Q) with respect to @ is defined in the 


previous development to be the matrix A. Because of the dimensions of the matrix A, 


the requirements for storage of the various matrices involved would exceed the 

capabilities of most computers considered for field usage. The structure of the 

matrix allows partitioning so that the sizes of the matrices involved and the 

computer storage requirements are both reduced. The niet A can be written as 
C 


A=3 


where the matrix C contains the partials of G(Q) with respect to the ground station 


parameters and B contains the partials with respect to the aircraft locations. The 


vector 6Q may then be partitioned into 


sq = ge 


od 


SEE 


with 6p the vector of station parameter Corrections and é6q the vector of aircraft 


location corrections. It is further assumed that the covariance matrix Ww; is of 


the form 


where wv and v are the covariance matrices of the apriori errors in the ground 
station and aircraft parameters respectively. This implies that there is no cor- 
relation between the apriori errors in ground station parameters and aircraft 
parameters which is a reasonable assumption. With these definitions the basic 
equation to be solved can be restated in terms of the separate matrices and vectors 
relating to the ground and aircraft portions of the problem. 

The original equation to be solved for 6Q is given by: 


(A7W A + W)8@ = ATW, 8R + W780 


so that with 
Ww 0 ° 
Cc 6 Pp ° 6p 
A =f . Q = [$2] ° W = ° 6Q = | 
Bj? ty 2 0 Wt}? e 
lI : : 
the equation becomes 
cw,C +H CWB 8 cw Wp 
1 Pp 1 P om sp +| P P 
Le Hu , ° 
BWC Bw,B+W | |6q) |B'w, W864 


The matrix Be is structured as follows 


T 
2 


vy 
-B 


% 


} 0 


a 


3(K-1) * J(K-1) J>3 


ag = Pygdiyy toe Gye Py, 
a 2 
Bie bps = Pa) ge ees MGp = Pp) E ae 


(agg - P3y)/Pyg vee» (435 - P3)/ry, | 3X7 


3 
R 2\1/2 
a > fag, EE) 
n=l] 


th 


Vz is the nth coordinate of the k” A/C location 


p.. as the n™" coordinate of the j*" ground station. 


nj 
Expanding the matrix equations into two equations gives, 
1) (c'w,c + wep + cl, BSq = C"W,6R + WP = 82, 
2) BW, Cop + (BW,B + W)8q = BTW, SR + War éZ, 
Solving the second equation for 6g gives 
6q = (BW ,B + " (82, - BTW, Cép). 


Substituting for 6q in equation 1) gives 


A 


" a] 
62 = (BW,B + Wy 82, 
Let the matrix V be defined by 
es ° oa 
v [¥, + c*y,c CTW, B (Bw, + Ww) B w,c| P 


Then 


6p = V7(82, - C7W,BE2,). 


Substituting back to solve for éq gives 


A -1,7 
6q = 62, - (B'W,B + WB Wy Cp 
The matrix Wy is diagonal with entries 1/0” so that the constant value 0? can 
be moved into the constraint matrices Wand W_ thus: 


T -1 = 
(c"w,c + W)6p + cly,B(B'W,B + W (62 - B'W,C6R) = 62, 


Z -1 # 7 A 
[w, + c"w,c = oT, BB WB+ ma) BWC |6p = 62, cTy,Bb2 


The vector ép is made up of 


“a 


Py ~ Py 
P21 ~ P21 
P31 ~ P3) 
where & is the apriori estimate of the 

parameter * being estimated from the data. 
Pie Sy 
ag = Pap 
8p = Psy - Psy 
of ~ 20; 
1g ~ 44 
2g ~ 224 


The vector 62, is given by 


62) = o*W ép + cl sr, 


P 
“C5 eo oe “Cy 6r 
Ts I 
« ow bp + -75 o ¢ @ t e 
-r o ¢ 6 Sh Sr, 
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The vector 62. is given by 


b2y.-% -B,or, + owe 8g 3x be, 


a 
825. = (BB, +07 We) M8209 Small ¢ 


27 
“C5 eo 8 @ Cy -Bo 825 
2 ; I ‘ 0 5 
BZ, = - : : 
ioe 2 e e e 0 
wed -By 8255 
ECB 82, 3a 
= ~EB 6255 J 
grip .éZ,.| J ] 
Tite 
$75 .82,..1 9 2 
tt 2t 
Then, 
62, = 62, - CBZ, = 0 “1 8 + [-2Cibr, - 3C;B.82,, | 
Sér. + SB Bn, : 
t 
-rr én, - cre. Sic j 
én. - rrlp 82 
t VIE® ACU Oy 
= oW bp + [-n0" (er + B.§Z,.) 
Pp t z t 20° t 
E( Sr ,+2B 625.) 
-59, (6r, . B82.) 


= are (8 rit B,8Z0,) 
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2 ° 
=o W dp + 37 
p?P 
J 
J 
J 
62, = br, + B.OZ,, Fx 1 


The vector 62) - c"Béz., is required along with 6z 


2° 


2. COMPUTATION PROCEDURE 
2 


1. Initialize a vector 6J x 1 witho W6,, (82,) 
2. Initialize a vector 3(k-1) x 1 with oH 8 (624) 
2 


3. Initialize a matrix to be diagonal o Wp to store the coefficient matrix V 
4. Start a "do" loop on zi = 2, k 


4.1 form matrix B, and save (J x WN) 


4.2 form matrix B, (W x J) : 
4.3. product 5B, (N x W) 


4.4 add o*w" to get BoB. + ofw® 
q GE 


4.5 invert (BB, + owe) and save 


Aa 


> 
(o,) 
ion) 
N 
TT 


825, - Bey. 
t tae 


2t 
a S Zt. 
4.7 82), = (BB, + OW) 829, 
4.8 62, = Or. + B 52, 
An A An i B. oF. T 
4.9 82, = &2, + [-Cf625, 625, -T)62,, -7,625] 
he 2,,t,-1 
4.10 product B, (BB, +0 i (J x N) 
T ted oe 
4.11 product B (BB, +o W) B. (J x J) 


eae axe 
4.12 form U, = I - B,(B;B, + o°W7) Be x J) 


4.13 form ru, * 
a & 4 

4.14 form ErU, 
4.15 form -£7.U. * 

tt 
4.16 form +£7U.?. a 

t+~v-I = 

e e * 

4.17 form rPU,?. { 
4.18 form Er,U,P - (J x J) ey 
4.19 form C7U, * (37 x J) a 
4.20 form -ECqU, (BI x J) 


A-10 { j 
| ; 


I 
I 
i 
i 
J 
i 
| 
q 
i 
4 
Al 
fl 
fl 
j 
A 
A 
i 
: 
i 


4.21 form ICyU 2. (37 x J) 

4,22 form Eee, * (Br x J) 

4.23 form IC1U.C. * (BI x3y) 
tee t 

4.24 form ClU.c. + , (BI x37) 
wee 


* indicates that individual matrices are formed by adding to proper locations in 
matrix V initialized above. 
4.30 end of do loop on iz 


1 


5. Invert matrix /V~ (67 x 67) 


6. correction vector $s -yl éz, (67 x 1) 
7. start do loop on i= 2, K 

7.1 form B; and save 

7.2 form B, 

7.3 product BB, 


7.4 sum BoB. + 0” 
Cae 


ye 
q ° 
7.8 invert (B°B. + 0“w*)~1 
tt q 
7.6 form vector [-c, ae r.| é,| (J x 1) 
7.7. product BY x vector + new vector 3 x 1 
7.8 product (BB, + ow)” x new vector 3 x 1 
7a 6q, = 625, - last vector Ks oa | 
7.10 end do loop on 7 
8. form new estimate of station vector = P - 6 
9. form new estimate of aircraft vector Q=@Q- 6 


q 
3. PROVISION FOR LOST DATA 

In the LRPDS problem it is possible that during the course of the flight 
terrain or aircraft maneuvers will shadow transmission between the aircraft and 
any of the ground stations. This will result in loss of range change measure- 
ments during times when the ground station is shadowed. Since the data loss 


may occur at different times for different stations it is necessary to make 


provision for an arbitrary loss of data for each station at each aircraft 
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sampling time. The method of handling the loss of data mst satisfy three 


criteria: 1) The method should allow the maximum utilization of the data collected 
so as to maintain the quality of positioning. 2) The method should not destroy 

the highly structured form of matrices involved so that full advantage of the 
structure can be taken in the computational procedure. 3) It is also necessary 

to retain the ability to provide a figure of merit for the resulting station 
locations determined by the systems which reflects the effects of the lost data. 

The approach taken to provide for lost data can best be visualized by a 
simplified example. Suppose that a linear least squares estimate is required of 
the n dimensional vector X from m measurements which constitute the vector Z. 
The observation equations constitute the rows of a matrix # where dimensions are 
mn, The matrix equation relating X and Z is thus 

HX = Z, 

The solution vector x then satisfies the normal equations 

HX = H'2. 

If in the measurement procedure some 2, (or set of 21) is lost or considered 
unacceptable there are two choices of how to handle the data reduction problem. 
First, the row dimension of the matrix H may be appropriately decreased and 
the measurements deleted from the vector Z. Second, a zero row vector may be 
inserted into H at each location of missing data and a zero inserted in Z 
corresponding to the missing data points. 

The second approach has the virtue of not disrupting any of the structure 
of H so that a fixed data reduction procedure will still accomlish the desired 
result. In order to accomplish this deletion of rows and data points in an « 
organized fashion define an m x m diagonal matrix S which has 1 or 0 entries as 
the particular measurement is present or missing. Then the relationship between 
X and Z becomes 


SHX = SZ, 


The solution X then satisfies the equation 
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I is J x J identity matrix 


T, is J«x J diagonal matrix with entries t, 


z is J x J diagonal with entries tape 
The vectors required are given by: 
. 2, 3 
624 csr +0 W5°P 


ent) 2 ie 
62 = BOR + ow bq 


fo) 
N 
I} 


T 2, \-1 
27 (BB+ OW) 62, 
bp = v-4(62, - c"x62,) 
1 2 
éq = 624 - (3"B + on) B Cép a 


The vector 6R may be partitioned into K-1 vectors each J in size, 


2 ory, 
6r= | ° ong, 
2 or; = : 
br, J(K-1) x1 y 
Srig| ox 
The vector 6q = 8419 ° 
8492 
8432 
Sa4x 
842x 


S4sx] 3(K-1) x 1 


s 


a 
ee 


vs | oy, +c [z_- B (BB + oy) ta | c| 
Sp = v1 (62, . c"86z.) 
and 
Sp = 2, - (BB + ofW)" B'Cép , 
with 
62, = Csr + on, ép 
62, = BSR + ow bq 


pede 22 mil 
ee sow.) 68, 


The matrix V which must be constructed is defined as 
vec [z oP a ow) ta" C+ on, 


Since B is block diagonal and on, is diagonal 


awe: (272, ae 


Ps Bb Hs ow) BY = 


U 


~~ 2.206). oF 
where U, = I - B, (B; + OW) E. . 


: Then the symmetric matrix V is given by 
: K K K 
ahs + ICU,C, - EC U, Sclu.r. scty.. 
t 4 t- tt tert a Se 
] 2 2 
i t=2 
K RO 
as BU.+ o°W -IU.P, -ru.t, ; 
t tt ey 3 
| 2 P 
EF. i a g 
te J ET.U.T. +0°W || OEPULT 
tTUL°t Pp tt 
3 
i Where we is the diagonal matrix of station location constraints Sry + owe 4 
Li : 
we is the diagonal initial range constraints ‘ 3 
| ee { 
Ww, is the freq. offset constraint ; 
| , ? , 
WP is the freq. drift constraint 
7 | The matrix C’ is given by } 
! > 
: 
7 ee | 
4 
ip ; ie” oh JE 4 
it Ce 
’ 1, “Ts sad Ty 2 
A ot te | 
i “Gy - Pape 
” 
, “Gay - Parry : 
“Gx - Psy 
| -@; ag Pad og 
ae ) 
- @az - P22)/r2; 


i] 
~-_ 
or 
we 
(] 
as 
ta 
Lo) 
= 
Ry 
° 
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ws? sux = H's"s2. 

Because S is diagonal and contains only ones and zeros 
sg=s 

so that, 
H’sHx = H' sz 


From this simplified example it is evident that the way to include provision 


for lost data points is to carry the effects of the matrix S through the entire 
data reduction procedure and delete those entries which would be multiplied by 
the zero entries in S. 

To relate the method of handling lost data directly to the LRPDS problem 
the development of the equations presented previously can be modified by 
replacing the matrices C and B with matrices SC and SB. The 6F vector is replaced 
with a vector SéR. The equations to be solved within the iterative portion of 


the solution method may then be restated as: - 


c" sw, 8¢ + C sw. SB 8 c'sw,s j 
Pp 1 Pp 1 ae PP 
. . 
B° sw, SC psw,sB + WJ 8 Bt SWS W 84 


Where all of the symbols are as previously defined. 

The primary influence of the inclusion of the lost data matrix S appears in 
the expression for the matrix V which is the coefficient matrix in the set of 
normal equations for the é5 vector. 


The matrix V now becomes 


K 
2 
o Ws + 26,02. “2040, oe fe 56,0. ost UR, 
t=2 
ow a sy. -br. -rU.7. “EU Rs 
p . ae 2 Lt D.% . 
V 2s 
G2? +.5ry.7. uk bt, orl R, 
2 ome e a. ie 
ory. 
6 WP +I? uh UR, 
of? + ER,U R, 


em ee ee — en es #26 fo 


Where U, = 5S, (r= BS, (B,5,B, + owt) $B .)8, 
5; is the eu J x J block (along the diagonal) of S, the lost data matrix. 

Another impact lost data has on the overall reduction problem is that if the 
number of ground stations visible from each aircraft location falls below some 
minimum value the solution for that aircraft location will either be impossible or 
so weak as to not be useful. If this happens then that aircraft location (or set 
of locations) must be excluded from the solution. 

The minimum number of visible ground stations allowable from each aircraft 
location is a function of how many coordinates of the aircraft position are being 
estimated at once. From geometric consideration, it is evident that for two 
coordinates being estimated (x, y) then there must be observation available from 
three or more non-colinear ground stations to resolve the possible ambiguity 
in the solution. Ina sintiar way, if three coordinates (x, y, z) are being 
estimated, then observation from four or more non-coplanar stations must be 
available. Although the use of four observations does not guarantee that the 
stations observed are not coplanar, it is certain that if only three stations 
are used, they will be coplanar (since three points uniquely define a plane). 
Without these constraints on the number of stations visible it is possible for the 
iterative solution to drive any aircraft point not having at least the minimum to 
an incorrect solution point corresponding to the wrong member of an ambiguous 
pair of solutions. 

4, REQUIREMENTS FOR SELECTABLE ESTIMATION OF PARAMETERS 

The formulation of the data reduction procedure up to this point has assumed 
that all of the parameters mentioned[i.e. 37 station locations, J initial ranges, 
J frequency offsets, J frequency drifts, and 3 (k-1) aircraft location parameters] 
are required to be estimated. Although this is true in the final stages of the 
adjustment, it is not true in the first few iterations to bring the parameters 
of interest into the range for final adjustment. The total data reduction program 
must therefore have provisions for selecting those parameters which are to be 


estimated at each stage of the data reduction procedure. 
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It is necessary to restrict the estimation to aircraft x, y only for the 
first few iterations in order to bring the relatively unknown flight path parameters 
close enough to the true values to perform the full adjustment on all parameters. 
This means that for the aircraft portion of the problem only two position parameters 
are to be estimated rather than three. In a similar way, the ground station para- 
meters to be estimated will vary depending on the degree of certainty of the 
solution. It, therefore, becomes necessary to provide for a flexible selection 
of the parameters to be estimated for both the aircraft locations and the ground 
station parameters. 
_— AIRCRAFT PARAMETER SELECTION 
For the aircraft locations there are only three choices to be considered: 
1) Estimate aircraft x, y and z 
2) Estimate aircraft x, y 
3) Do not estimate aircraft locations. 
— * GROUND STATION PARAMETER SELECTION 
For the ground station parameters the number of choices is more complex. 
aout Station Position Coordinates 
1. Estimate all ground station position parameters - x, y, Z 
2. Estimate ground stations x, y only 
3. Do not estimate ground station positions 
Ground Station Initial Range - ao. 
The estimation of the Qos term for each ground station is required for 


all circumstances. 


Ground Station Frequency Offset - 245 

The ground station frequency offset, ay; may or may not be required 
depending on the level of accuracy required at a given stage of the data 
reduction. 


Ground Station Frequency Drift Rate - aos 


The frequency drift rate, ayss must be selectable. 
A-18 


5. DATA REDUCTION FLOW DIAGRAM 
The overall requirements for the data reduction and the equations to be 
mechanized have been established by the preceding discussion. The outline of the 
data reduction procedure may now be formated to include all of the provisions 
required. 
In order to derive the flow diagram for the computation the following 
notation is used: 
J - number of ground stations being processed 
NG - number of ground station position coordinates to be estimated (2 or 3) 
(et Yi OF XK V4 2) 
IGC - control parameter (1 or 0) 
IGC = 1 enables the estimation of the ground station positions 
IGC = 0 disables the estimation of ground station position 
Pag 3 m= 1, NG; j «1, J) 
IAl - control parameter (1 or 0) 
IAl = 1 enables the estimation of ground station frequency offset 
(33571, J) 
IA2 - control parameter : 
IA2 = 1 enables the estimation of ground station frequency drift 
rate. (a955 j= 1, J) 
K - total number of aircraft locations used 
NA - number of aircraft position coordinates estimated (2 or 3) (x, y or 
Xy ¥» 2) 
IAC - control parameter (1 or 0) 
IAC = 1 enables the estimation of the aircraft location 
IK(i) - control parameter 
IK(i) = 1 enables the estimation of the jth aircraft location 
IK(i) is set to zero if the number of ground stations visible 


fails to satisfy the necessary criterion. 
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NSTA - minimum number of stations required to be visible from each aircraft 
location 
NSTA = NA + 1 
IDATA (j, i) - Control variable (1 or 0) which defines the diagonal elements 
of the lost data matrix S. The existence of a 1 at location 
j, i of the IDATA (j, i) control matrix denotes that the 
j= station was visible from the ith aircraft location and 
conversely if a zero exists. 
IACS, IGCS - Control parameters (1 or 0) 
IACS = 1 or IGCS = 1 denotes that one or more aircraft locations 
or ground station parameters have not satisfied the convergence 
criterion and hence another iteration of the problem is required. 
When all parameters have satisfied the convergence criteria then 
IACS and IGCS will be zero. 
Jl, J3, J4, JS - indices of the location of various parts of the matrix V. 
The symmetric matrix V is organized in block form as shown 
below 


V(s1, J1) 


V(J1, J3) 


V(J1, J4) 


» V1, JS) 


NG*J x NGS] NGS x J | NGI x Jo] NGS x J 


V(J3, J3) | V(J3, J4)| V(J3, JS) 


J x J J XS 


V(J4, J5) 


V(J4, J4) 


JxJ. Jx J: 
V(J5, JS) 


“JIxJ 


The indices shown refer to the first TO’ and column (upper, left hand corner) of 


each sub-matrix. 


In a similar way the vector 6Z 


62, (31) 
82, (33) 
62, (34) 
62, (J5) 


1 is denoted by 
NGeJ x<1 

Jx1 

Jxil 


Jxi1 


to distinguish the various component parts of the vector. 


In general 


J1 - defines the location containing elements of the matrix (or vector) 


related to ground station position. (Prj) 


J3 
J4 


defines the location related to initial range estimation. (a9;) 


defines the location related to frequency offset. (a,,) 


JS - defines the location related to frequency drift rate (a5) 


Si; $2; 


S3, S4, S6 - Refer to temporary matrix storage locations required 


in the evaluation of the matrix V. 


DRP - is a temporary vector storage location 


ITIME - 


control variable (1 or 0) 


ates 


The use of ITIME allows a portion of the program to be used for 
two different purposes. When ITIME = 0 the path through the 
computation procedure yields the matrix V and the vectors 62, and 
éz, which are used to compute 6p. ITIME is set to 1 and a second 


pass made. through the program to compute the aircraft location cor- 


rection 6q. 


Notation: The notation AtA + B means that the new value of the variable A is 


determined by:iadding B to the eld value of A. In most cases this 


operation could be written as A; = A; _ , + B;, but in order to 


1 


reduce the notational requirements the form listed above was used. 
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APPENDIX B 


COORDINATE CONVERSION CALCULATIONS 


1.0 INTRODUCTION 

The requirements for the LRPDS program state that positional 
data input to and output from the computer terminal be in either 
Universal Transverse Mercator (UTM) or Geodetic Latitude and Longi- 
tude Coordinates (LAT-LONG). The nature of the data reduction 
problem is such that it is highly desirable to utilize a simple 
local X, Y, Z coordinate system thus avoiding many trig function 
calculations. The most direct way to approach this coordinate 
conversion problem is to convert all coordinates to lat-long and 
then from lat-long to the coordinate system desired. This method 
also allows for easy solution to the adjacent zone problem in 


overlap regions since conversion to Easting and Northing in each 


are required; namely, LAT-LONG to UTM, UTM to LAT-LONG, LAT-LONG 
to LOCAL, LOCAL to LAT-LONG. 


2.0 THE SPHEROID 
The accepted shape of the earth for mapping purposes is a 
spheroid generated by the revolution of an ellipse as shown in 


Figure 1. 


en 


i zone can be made from the lat-long position. Thus, four routines 


Figure 1. Meridian Ellipse of the Spheroid 


In Figure 1, AP = N is the normal to the meridian ellipse at 
P(r, z). The semi-major axis is a and the geodetic latitude angle 


is $. The equation for the ellipse is 


r2(1 - e7) + 2? = a? (1 ~ €?) q) 
where 
2 
e-  e e) (2) 
a 


Sp 


pate 


and € is the eccentricity, b is the semi minor axis of the ellipse. 
The spheroid specifications are usually given as t..e dimension of 


a and the inverse flattening 1/f. This factor is related to the 


other parameters as 


eh x 


f a-b 


b = acl - f) (4) 


e = #0. 3) (5) 


The list of spheroids in common usage are listed in Table l. 


Since the slope of the normal at P is the negative reciprocal 


of the slope of the tangent. 


dr 


tan? = = >a 


z= r(1- €”) tand = N(1 - €7) sing 
Substituting the value of z in terms of r we have 
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The radius of curvature of the meridian ellipse is given by 


a(l - €7) 


2 


QQ. - e? sin? $)3/2 


(10) 


In Figure 2 the ellipse has been revolved about its minor axis 


through an angle A with the point P moving to p!. 


coordinate system, we have 
y =r sin A = Nos ¢ sina 
x = Yr cos A = N cos ¢ cos A 
Z=N (1 - «2) sin ¢ 


x2 fs y” ze - 
ee ee aes 
a 


Zz 


In the X, Y, Z 


(11) 
(12) 


(13) 


(14) 


Figure 2. Generation of Spheroid 


3.0 THE UNIVERSAL TRANSVERSE MERCATOR GRID AND PROJECTION (UTM) 

The transverse Mercator projection is a conformal projection 
in that the angles measured on the projection or computed from the 
grid coordinates closely approximate their true values as well as 
the local scale factor at any point is the same in all directions. 
If we consider that a Mercator projection is similar to a projection 
upon a cylinder encasing the earth, the orientation of the cylinder 
axis for a transverse projection would be in the equatorial plane 
with the cylinder tangent to the earth at the central meridian of 
the map projection. 

There are 60 transverse Mercator zones each 6° in width extending 
north to 84° and south to 80°. Zone number one lies between 180° 
and 174° west longitude. The zones are numbered consecutively 
from west to east. A 50 mile overlap provision is provided for at 
each zone boundry. | 

A scale factor of 0.9996 introduced along the central meridian 
of each zone gives the effect of a secant condition in the geometric 
sense. 

The zone contains a metric grid superimposed upon it with a 
500,000 meter false Easting along the central meridian, a zero meter 
Northing at the equator for the Northern Hemisphere and a 10,000,000 
meter false Northing at the equator for the Southern Hemisphere. 

The ellipsoid on which the projection and grid are based 
depends upon the area involved. The Clarke 1866 is used in North 


America. Other areas are as listed in Table l. 
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3.1 UTM GRID FROM LAT-LONG COORDINATES 


Figure 3. The Linear Element of the Spheroid 
is Obtained from A Differential 
Right Triangle 


The transverse Mercator projection is a conformal mapping of the 
spheroid onto a plane such that a rhumb line intersects each 


meridian at a constant angle. From Figure 3, we have 


tana - —Ncos ¢ da (1s) 
rd¢ 
or 
R 
da tan a@ + —y- sec gd d¢ (16) 


Setting @ equal to a constant the integral curves are 


A- XA}. = Ttana (17) 


where T , the isometric latitude, is given by 


t= if —~ sec @ dod (18) 


°o 


now in mapping the z plane onto the w plane, 
w=x+iy = £(z) = f(A + iT) (19) 


The requirement for a transverse Mercator projection is that the 
scale shall be true along the central map meridian. Hence when 


A’ = 0 we must have x = 0 and from (18) and (19) we must then have 


ly = f(ir) =i S¢ (20) 


where S% is the arc length along the elliptic meridian of the 
spheroid from the Equator to latitude $¢. 


3% * f° Rd¢ (21) 


and from (18) we have 
Rd? = Ncos ¢ddtT (22) 


So that 


i 
Sy = * N cos @¢ dr = f(T) (23) 


= 


If we expand x + iy = f(\ + it) about the point z = it we obtain 


' 2 ‘ 
x + iy = £(A + iT) = £717) +ALL (it) + Ee f "(ir) 


fit (it) + erececvece 
ramet) + eeer0ee#e? a 


now from (20}and (23)it is seen that 
F(i7)-= 184:= 4f(1) (25) 


differentiating this expression and substituting into (24) we have 


x + 1 = if(T ’ 2 ef 2 ete 
y De geen REDS Si Set oe WON rea ot (7) 


vili 


if (Tv) eae tte 


equating real and imaginary parts 


»» 


5! 


’ 2 eee 
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2 4 
y= f(r) = = f'' (1) + AL 
2! 4! 


+ 


Using the following derivatives and trig identities: 


N = (N - R) tand 


R= 3 (N - R) tan¢d 
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N 2(N — R) 

| acs ae meal 

qd? _ UN 

ar i‘ “R_ Eop @ 

(N cos ¢) | = -R sin¢g 

(N sin ¢) oe Sec¢d (N-R sin? @) = (R cos ¢)/(1 - e*) 

2 sin n@d cos @ = sin (n + 1)? + sin (n - 1)¢ 

2 cos n¢@ cos ¢ = cos (n + 1)¢% + cos (n - 1)¢ 

2 cos n @ sin ¢ = sin (n + 1) ¢@ - sin (n - 1) $ 

2 sin n@? sin ? = cos (n - 1) @ - cos (n + 1) ¢ (28) 
Setting N/R = o0 we have, 

f' (tT) = N cos od 

a (tT) = cnt cos @ - N sin ¢) ch a = cae sin 2¢ 

“dT <7 
(7) = we [are y cos ¢@ + (7 + 1) cos 3¢ 
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wit = 2 


2(1 - 2@ + 1307 - 40°) cos ¢ 


+ (-3 + 20- 307 + 440°) cos 3 > 
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+ (1 + 20- To + 2807) cos 5 ¢ 


pn eee Swen oer oom ee, 
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-~ 170807 + 35920”) cos 5¢ | 


Substituting (29) into (27) we have 
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o- N-14 8 cos” (31) 
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We may rewrite (30) as 
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Arranging (35) and (36) to facilitate computation and setting 
AK=) = Xo = AX as the longitude difference from the central 


meridian and dropping out small terms we have 
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The first three terms of (37) correspond to the functions IV, V 

and Be given in the Army Tech Man TM 5-241-8. The first four terms 
correspond to the functions I, II, III, and Ag of the same manual. . 
However, the functions Ag and Be of the manual are truncated ver- 3 


sions of the terms shown above. 
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3.1.1 Meridian Arc Length from Equator 
The meridian arc length § is defined in equation (21). 


Substituting the value of R from (10) into this equation, we have 


) 3/2 
Sy - JS ails <°) @ =<" wie) dt (39) 
This may be written as 
Sy- a(l- 7) I (, - €7, €) (40) 


where J7(¢, - Ca €) is Lengendre's Normal Elliptic Integral of 
the third kind. Now 


(o, ~ €7, €) = Saree (6, €) - 


(1 - €? gin? 6)7)/2 ee? sing cos $) (41) 


where E (¢, €) is the incomplete Elliptic Integral of the second 
kind. 


A series expansion for So is given in Jordan's Handbook of 
geodesy, Vol. III, first half. This expansion using all terms 
involving e« through the 10th power is 
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3.1.2 Central Scale Factor 

Since the local scale factor increases as a function of the 
distance from the central meridian, the UTM system has adopted a 
central scale factor of 0.9996 to equalize the scale errors across 
the zone width. This scale factor is the one used for this program. 
Two other scale factors are in use in other areas of the world, 
such as .9999 used in Canada and 1.0000 used in many European 
countries. The scale factor is most easily accommodated into the 
calculation by modifying the semi-major axis by the scale factor 


as 
' 17 
a= k,8 (44) i ; 


where a’ is the value to be used in calculating the UTM coordinates 
and Ky is the central scale factor. The units of a used in these 
equations is in meters instead of the usual units of minutes on 


the Equator. 


3.2 LAT-LONG COORDINATES FROM UTM GRID 
Rewriting (19) as 


AX+ it = F (x + iy) (45) 


Now when x = 0, A= 0 and then F(iy) = i from (22) and (23) 
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Expanding F(X + iy) about the point iy r 
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From F(iy) = it , we have F' (iy) = r', F" (iy) = -it", 


° 
if 


Fiii (iy) = -riii, FLY (iy) = iriv, eee pn’ = i(-i)® a t 
We have . 
2 ; 3 
A+ it = at) + a7," - Se it, Og ie 
n t - 
eece = i(-i) "1" + Oi OTs (50) 
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where the subscript one refers to the "footpoint'"' latitude of a 


given point of the projection as shown in Figure 4, =) 


| 
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Figure 4. Footpoint Latitude %, of Point P 


from (47), (48) and (28) 
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- 1536t1n* ~ ssoan® + 27456¢7n® 4 744t4n® 


10 


114727)° 4 53952t2n! ~ 7872+4n8 


467277 4 30528t-nl2 ~ 2ag60t4n!® 


\ = 23040t4n12 
N 


‘“ 


| 
7 
a es 


Placing these values into (51) we have 


ee 


x 1 x - 2 2 
[7 - (¥) SE eats tae 


1 
: 2 2 4 2. 2 ) 
(+ wollte \ 
| 4 6 2, 4 2,6. 
Su 24t, - 47, + 4t, 1 + 24t, " } 
\ 
\ ~ 
i 61 + 662t.7 + 1320t.4 + 720t.& —\ 
\ 4 os 1 1 \ 
Aipethe ! + 1077.7 4 437.4 + 440t,20_? 
: 1 1 lg | 
cos . / 
a 2, 4 a2 ol 
| _ \+ 977° ~ 234t,°"* + 336t, 9, | 
i 4 
Ny + 1ss7,2 _ 7724, % ~192t.47,.4 | 


1 


4, 8 
1536t, i iis 1632t 


11 


+ 

S&S 

3 

_ 

°o 
ial ace 


H 
4, 10 - 
1920t, 7) - ; (60) 


eS. 


7 
In the coefficient of (+=) , all terms containing 7, may be 
1 


deleted without seriously effecting accuracy. 


; 61 + 46 n° 


2 
F 3 6 + 48 t) 1 
zr 
2,4 
- 36t ) un 


2 
- 192t, nN 


if 
Nan 1o920t, 4 


6 


i 
+ 5040t, 


6 
- 29277, 

8 
- 88087, 


- 114727, * 


- 7872t, 97 


4 
- 24960t, 


Since Avr is the difference of the true isometric latitudes in 


order to obtain the geodetic latitude difference from (61) we 


expand Ad into a series in AT as 


[ 1385 + 7266t, 


+ 1sot,” - 


2 4 6 


+ 120t, + 100%, 


2,, 6 8 
- 96t, 1 


+ 887, 


8 


aoe 17317,” 


2 


2 
+ 4416t, an 


4 


& 1830t,°7 4 2 


4 
7 + 2688t, 1 


256 4,4 
+ 5052t, Ls - 1536t, ua 


4 6 


2., 8 
+ 744t, un 


+ 27456t, 1 


0 2 


+ 53952t,°7 10 


1 


8 
1 


2 . 30528t.27 


1 
- 4672", i %4 


10 4 


- 23040t 12 


n 


q eas 


1 


CL ee eS ee 


a 


i 
I 
if 
i 
B 
l 
l 
I 
E 
b 
; 
f 
I 
| 
| 
| 
| 


2 
d¢, 2 d“¢ 
Ad $ $, AT “oF, # 7 are 
1 
At"a"d, 
AOS eee es eeoeeee (62) 
n! at 


Now from (46) 


sh “ ws cos $@ = (14+ 7 2, cos ¢ 
Ri 1 1 1 


T 
a7) 


and by successive differentiation 


See es Ae 7,7) G« an yt, cos” $, 


—- = (14 0,7) [ 40)? + 200 + 02700 | 9t,?) 
2.2 2 3 
+ 3(1 + 7,°)" (5t,” - 1)| cos $, 
a = ty (8) | -st,? + 41 + 9,7) (21t? = 4) 


+ 4(1 4 7,7)? (17 = 45t,7) + 30 4 Nw 


2 4 
(35 t,~ - 19) | cos “ 9, (63) 
ano, 
Substituting the value of Ar from (61) and the values of A 
at. 
1 
from (63) into (62) and deleting terms in uy in the coefficient 


8 N 
of (+) as well as taking note that 1 + :° = a we have 


1 


1 4 2 2 
$ % -ory * tg x (5 + 3t, + 
191 24R,N, 
2. 2 
- 41,7 - 9n,7t,7) 
2 2 4 2 2 
61 + 90t,” + 46,” + 45t,* - 252t,°n, 
t 
_ —— 7 = 3n,4 + 1007, ° = e6t,°n,4 ee got, 4n,? 
720R,N, 
8 4 4 2.6 2 8 
+ 881, + 225t, Y + 84t, " - 192t, " ! 
ty 8 2 4 6 
+ ———__7 x {1385 + 3633t,” + 4095 t,7 + 1575t, (64) 
40320R,N, 


Equation (60) with the terms containing at in the coefficient of 


(-)" deleted, along with equation (64) are used to determine 
the , from the given X, Y location. The AX is added to the 
value of the central meridian to determine A. The X and Y values 
used in the equations have the false Easting and false Northing 


removed. 


3.3 COMPUTER CALCULATION OF UTM FROM LAT-LONG 
This routine required double precision floating point arith- 
metic on a Sigma V computer to maintain sufficient accuracy. An 


outline of the program is as follows. 


a. Input 
Semi major axis in meters =a 
Central scale factor = kg (.9996 for UTM) 
Inverse flattening factor = 1/f 
UTM zone number 
Latitude = dradians 


Longitude = A radians 


b. Step l 
Calculate 
a ke & 
f = (1/(1/f)) 
e2 = f (2 - f) 
ha = (~183 + 6 (zone No) )/57.29.... 


AA= A- Ay 


2 1/2 


N= a'/(1 - €? sin? 4) 


c. Step 2 


Calculate S¢ from equation (43) 


d. Step 3 
Calculate x from equation (37) 


y from equation (38) 


Step 4 
Add 500,000 meters false Easting to X. Test Y and if 
negative, subtract Y from 10,000,600 meters false Northing 


(i.e. add magnitude of Y to 10° meters). 


3.4 COMPUTER CALCULATION OF LAT-LONG FROM UTM 


This routine required double precision floating point arithmetic 


on a Sigma V computer to maintain sufficient accuracy. An outline 


of the program is as follows. 


a. 


Input 

Semi major axis in meters = a 

Central scale factor = Kg (UTM, Ko = .9996) 
Inverse flattening factor = 1/f 


UTM zone number 


X = Easting containing false Easting (meters) 
Y = Northing containing false Northing (meters) 
Step l 


Remove false Easting from X 
Test and remove false Northing from Y 


Ao = (-183 + 6 (zone no))/57.29.... 


a = Kg a 
eo =F ee, = 
Step 2 


Iteriatively solve for footpoint 


Latitude b, by the following method 
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ce, Calculate Coefficients 
Ca through Ce of eg (43) | 


Aye 
iso | G, = AY/ (Ca (1-¢°)) 


Solve Foy AY = Sy. 
using oy (43) 


Step 3 
Solve for Ad from equation (60) 


Solve for $ from equation (64) 


Step 4 


Solve for A as Xo + A 


4.0 LOCAL COORDINATE CONVERSION 
4.1 LAT-LONG TO LOCAL COORDINATES ; 
The local coordinate system chosen for data reduction is a 
cartesian X, Y, Z coordinate system with the Y axis oriented north 
at the origin, X oriented to the east at the origin and Z vertical. 
The X-Y plane is normal to the spheroid at the reference or origin 


point. 


Setting up a geocentric coordinate system with the X, - Xg 


plane in the equatorial plane and X3 oriented north as shown in 


Figure 5, the coordinate system x, as Shown in Figure 6 is related 


to X, as follows: 
¢ t 
X, = X) cos A- Xo sina 


¢ 
Xo = X) sin A - 2 cos Xx (65) 


and 


1 
| 
; 
| 
| 
| 
i 
| 
{ 
I 
j 
i 


Xo = (0) (66) 


' 
X, = E sin o, +h sing 


where \ is the longitude angle, ¢ is the geodetic latitude, es 
is the geocentric latitude and h is the geodetic elevation 


of point M. 


B31 ; 


Yeo Mende 
a 


Figure 5. Geocentric Coordinate System 


Equation (65) also may be written as 
X, = Cos d (E cos %, +h cos >) 


X_ = Sin ) (E cos an + h cos ¢) (67) 


X3 = E sin ¢,+hsing 


5 


. 


' 
Figure 6. Rotated Xk Coordinate System 


Now on the surface of the spheroid, the intersection of the 


' 


xy Xo! plane with the spheroid is an ellipse expressed by 


(68) 


A lant BS, BO, 


mn a 


B-33 


So that 


ee OS OR ET, PURI HE 


2 een d E2 adn i) 
E Bors. 4 29) (69) 4 
a b \y 
; 
or | 
E = ig ag (70) 
(cos $,)/a + (sin $,)/b ; 
Since ¢ is the angle of the normal to the ellipse and 
dx") 
tan @¢d = - : (71) 
dx 3 | 
differentiating (68) with respect to x3' 
, ' s] 
xy dx j X 3 : 
ee eee ee = 
a dx 3 b 
So that 
az x! ar 
tan? = — ; i — tan d (73) 
b xX 4 b 
Since e2 = (a2 - b”) /a” 
tan @ = (1 - e7) tan¢ (74) r 
and 


a: ae eet ge — 


- 


4.1 LOCAL COORDINATE SYSTEM 
Let ae Xo" and oe be a coordinate system with its origin 


at reference position A as shown in Figure 7. 


, 
Xq 2%, hy 


x, aud x," 
directed we 
Paper. 


%, 


Figure 7. Local Coordinate System with Origin at A 


Then 
x,"| 0 1 0 x, - x4 
x2" | = '=j sin da 0 cos aN Xo (80) 
po cos dy 0 sin os Xo" - Xs 
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Se EE LIENS 


X) - Xz, 4 cos Ay sin ra 0 x, - X14 
t 

Xo = - sin An cos Ay (0) Xo - Xon (81) 
’ ) 

X, - Xo 4 (0) 0 1 X3 - Xgq 


x = AD (82) 
where 
*E “ in 
D = Xo - Xon (83) 
my” ROR i 
and a 
- sin Ay cos An 0 
A = - cos h, sin >, -sin bs Sin h, cos $, (84) 


L cos 6, cos hy cos $, sin \, sin ?s 


The local coordinates can then be determined by use of equation 
set (67) to solve for X and X, utilizing (74) and (75) to determine 


$y> tga» E and Ey. 
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\ 


Sara ae : 


A ee 


4.2 LOCAL COORDINATES TO LAT-LONG 

The inverse solution is not as straight forward as the LAT- 
LONG to LC in that the set of equations (67), (74), and (75) 
relating X to LAT-LONG-h are not linear. To solve this set, some 
sort of iterative approach is necessary. One method would be to 


rearrange (82) as 
x =A x + Xq (85) 


and then iteratively solve for LAT-LONG from x. This, however, 
introduces some error in the calculation in that errors in the 
inversion of A are propagated to the answer. A better approach is 
to estimate the geographic coordinates and solve for the LOCAL 
coordinates in a manner to minimize the square error between the 
input local coordinates and the computed local coordinates derived 
from an estimate of the geographic coordinates. Figure 8 gives a 


block diagram of the solution where 


Geographic coordinates are 
Wy = @, Wo = d, W3 =h (86) 


Local X Y Z are 


1 =X, Xo = JY, ~~ = 2 (87) 
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eres 


ORS: OTE 


ws estinated wW 


Figure 8. Inverse Solution Block Diagram 


In this solution, the block relating x" as a function of W is the 
LAT-LONG to local coordinate conversion routine outlined in 4.1. 
Matrix C relates the geographic error vector to the local coordinate 


error vector in a manner to minimize the square of the error between 


we A 
x" and x". The vector equations are 
T Pr, T " 
BB bw =B Sy (88) 
or 
a Tern =]. eT owe " 
by (B’B) B 5, Cc oy 
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From (82) and (83) 


a1) 
x = Ax - AX, 


So that 


B= -~AF 


where A is defined in (84) and 


Ox, 
OW, 


Ox, 
OW, 
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Ox, 


OW, 


OXo 


OW, 


Ox, 
OW, 


(89) 


(90) 


(91) 


(92) 


Reaves. 


Since 


t, =¢ (93) 


we may approximate (67) as 


| 


xX) = (E + h) cos @ cos A = (E + Wo) cos (W,) cos (Wo) 


5 (E + h) cos @ sin dA = (E + Wa) cos (W,) sin (Wo) (94) 
Xq_ * (E +h) sin @ = (E+ W3) sin (W,) 
Then 
Ox, 
aW, = - (E + Wg) cos (Wy) sin (W) 
Ox, 
aW, = - (E + Wg) sin (Wy) sin (W,) 
Ox, 
ow, = (E+ W3) cos (W,) 
Ox, 
BW, = - (E + Wg) cos (W,)) sin (Wo) 
OX5 
a, =  (E + Wy) cos (W,) cos (Wo) 
Ox, 2% 
OW, 
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cos (Wo) cos (W,) 


OW, 
OXo 
OW, = sin (Wo) cos (W,) 
OX, 
= sin (Wy) (95) 
dW, 


4.3 COMPUTER CALCULATION OF LOCAL COORDINATES FROM 
GEOGRAPHIC COORDINATES 


The computer program utilizes equation sets (67), (74), (75), 
(82), (83) and (84) to determine the local cartesian coordinates. 
Although some problems can occur in solving for points close to 
the coordinate origin where the major positional difference is 
altitude, utilization of double precision arithmetic on the Sigma V 
computer appeared to be quite adequate for this method of solution. 
This avoided the problem of requiring different routines for dif-~ 
ferent regions. 


The method of solution is as follows: 


Step 1 
Input 
a = Semi major axis in meters 
1/f = Inverse flattening factor 
R; = Reference Lat-Long in radians, altitude in meters 
Wy = Geographic coordinates of point to be converted in 


LAT, LONG in radians, altitude in meters 


Step 2 


Calculate 
f = ACESS) 
«2 = £(2 - f) 
b = a(l - f) 
E, = b/(1 - €? cos? $4) 
bea = tan? ((1 - 7) tan $,) 


Set up A array equation (84). 
Calculate 


X, as RX, from equation (67) 


Step 3 
Calculate 
E = b/(1 - €? cos? ¢$) 
by = tan4((1 - e”) tan 6) 


x as PX; from equation (67) 


4.4 COMPUTER CALCULATION OF GEOGRAPHIC COORDINATES FROM 
LOCAL COORDINATES F 


The program follows the method outlined in 4.1 utilizing the 
conversion method of 4.3 to convert estimated geographic coordinates 
to local coordinates. Double precision on the Sigma V is utilized 
in the arithmetic and appeared to be more than adequate in the 
solution. A test of all programs was run utilizing the UTM- 
GEOGRAPHIC routines as well as the LC-GEOGRAPHIC, This test in- 
volved setting a reference on a reference meridian and offsetting 
the point in UTM grid by increments of 10 meters in Northing and 
Easting and varying altitude from 0 to 500 meters in steps of 100 
meters. The UTM points were converted to Geographic to LOCAL to 
GEOGRAPHIC to UTM with an error of less than .0001 meters altitude, 
.001 meters in X and Y and 1.0 x 107° radians. 


The inverse solution method is as follows. 


Step 1 
Input 
a = semi major axis in meters 
1/f = inverse flattening factor 
Ry = reference geographic coordinates in @?, A, h 
x" = local cartesian coordinate value of the point. 
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Step 2 


Calculate 
f = 1/(17f) 
€? = £(2 - f) 
b = a(l - f) 
E, - b/(1 - €? cos” 4,) 
$e, = tan”) ((1 - €7) tan $4) 


Set up A array equation (84). 


Calculate 


Xq as RX, from equation (67). 


Step 3 


Choose initial we = R, 


Step 4 
Iteration procedure as in Figure 8 
Calculate 


b/ET =" aoa BY 


Leo | 
" 


$, = tan”) ((1 - €7) tan 


x as PX from equation (67) 


Step 5 
x" = A (x - X,) 
Ax" = x" = ¢" 
Test Ax 


If Ax < .01 stop iteration. 


Step 6 
Set up F array as in (92). 
Calculate B and B! as in (91). 
Calculate (B’B)~! pl = a 


Calculate new we as wo-A Ax" 


Return to Step 4. 
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APPENDIX C 
PROPAGATION CORRECTIONS 
1. INTRODUCTION 
The propagation correction procedure discussed in this appendix is a 
simplified version of that reported in the LRPDS Error Analysis dated 15 April 1971. 
The procedure is simplified to the degree that each of the 1800 possible range 


ibta Sigma A Eo i 


measurements made during a flight can be corrected individually if desired. 


i i Certain simplifying assumptions are discussed in the following paragraphs 
that allow the above conclusion to be made,. 
. The National Bureau of Standard at Boulder, Colorado recormends (reference 
.& 
. & 


one) the use of a bi-exponential model for the refractivity of the atmosphere, 
N(h) = D, eng (7m +W, en Wha (1) 
in which the first term is referred to as the "dry term", and the second is 
called the "wet term". In equation (1), h, is the surface altitude in kn 
above mean sea level at which the coefficients D A and W g are measured; the 


latter are determined from the Smith and Weintraub constants (reference two), 


is ns P 
| (2) 
: W. = 3.73 x 10” © 


8 


A me RRP A 
“ r 


where P is the total atmospheric pressure in mb, 7 is the absolute dry bulb 
temperature (°K) and e is the partial pressure of water vanor in mb. The 
partial pressure of water vapor may be obtained from.the psychrometric formula 
(reference two) 

e = e% - P(?-T) (471 + 0.665) x 10°° (3) 
where e; denotes the saturation vapor pressure in mb at the absolute wet-bulb 
tempterature, 7°. It is assumed in this discussion that the measurements 
required to make an accurate determination of D, and W, will be available. 


It is further assumed that reasonable estimates for the values of the dry 


and wet terms in the vicinity of the aircraft altitude can be made; from 


e*T 


pls alle ecard Batam ia 
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this information, the scale factors 6 and w in equation (1) can be determined 


as 
5 = PnlD,/0(h)] 
Th 
(4) 
and yp = 27lW,/w(h)] 
‘ m, 


Given that the atmospheric refraction is modeled as a bi-exponential, what 
are the residual propagation errors caused by modeling errors? This question is 
answered in reference three. In that study, both horizontal and vertical 
refractive gradients were assumed and the effects on the value of the range 
correction were determined; the horizontal gradients were approximately 60 N- 
units per 80 km. The NBS people at Boulder suggest that this represents the 
worst case Situation for typical locations around the globe. As many as five 
different vertical profiles were inserted into a computer program that evaluated 
the actual refractivity along each ray path hetween an aircraft at 6 km altitude 
and FO's located at various points on the ground. Assuming that this atmosphere 
could be represented in bi-exponential form, with the parameters in equation 
(1) known via two measurements, one made on the ground and one made at the air- 
craft altitude, the results showed that the ms variation in the range cor- 
rection due to horizontal gradients was about 1 part in 10°, 

However, the above assumes that a measurement would be available in the air- 
craft; no such measurement will be available in the field. Reference four showed 
that if the refractivity at the aircraft (or any other arbitrary) altitude is 
estimated rather than measured, then the loss in accuracy is about 1 part in 10° 
for every 10% error in the estimate. This error proves to be the dominant , 


residual error assuming that the surface measurement of refractivity is good to 


within 5% (typical radiosonde measurements are better than this). 
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OORRECTION PROCEDURE 


It has been shown that the atmospheric refractivity can be modeled, with suf- 
ficient accuracy, as a bi-exponential. A sufficiently simple procedure to obtain 
the actual range correction is next described. 

The effects of an inhomogeneous atmosphere were examined in reference three; 
there it was shown that inhomogeneities in the refractivity were averaged over 
the ray path. Specifically, large refractive gradients (including surface ducts) 
had very little effect on the range correction required; the range correction 
was sensitive to only the average index of refraction, . Therefore, a simplified 
procedure requires a means of determining 7. Assume the typical propagation 


geometry of Figure C-1. 


Figure C-1. Typical propagation geometry for which h, in the surface altitude 
above sea-level, a is the earth's radius, 6) is the geometric angle 


of elevation, and 6 is the angle between the ray P,P, and the normal 


to the radius a + h. 
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It is seen that the actual path length, *%,, between points P, and P2 is | 
given by ‘ 

R_ = car 

a 
where @ is the average propagation velocity and AT is the existing time delay. 
However, @ = e,/n where eC, is the propagation velocity in a vacumm and 7 is the 
average refraction index. Therefore 

e 
R,*oaT= 1k, 


n n 


where R, is the "measured" range. By definition, the range correction is the 
difference between the measured and actual ranges. 


AR=R -R 


m a 
m (fie B) Re (S) 
= (n-1) R i 
n 
By definition, the average index is given by* = 
P2 raat 
oe i) W(R)dR (6) i 
Ge P, - 
along the ray it is seen that dh = dR sin 6, so i) 
: h ee _ 
> al f° Bet a (7) { 
Ry sin 6 
h 2 


§ 


*It is assumed that the geometric path P,P, does not anpreciably differ from ji 
the actual path which is slightly curved due to refractivity. It has been 


shown that these two paths differ in length by only centimeters for typical 


applications; consequently, the expression in equation (6) is essentially 


correct. 


——E 


> 


where .n(h) = 1 + W(h) ° 10°; N(h) is given by eqn (1) and is measured in 


N-units. Breaking the two integrals into two integrals, 


h -6 h 

ee a 10 ‘a N(h)dh 

ner dh + » 

Lal f aa R, f sino (8) 
h h 
8 8 


it is seen that the first integral is first equal to Ro Therefore, eq (8) 


becomes 
h 
-6 a 
R, Sin 6 


It has been found that the average refractivity given by eqn (9) is remarkably 
insensitive to changes in 6. For instance, assuming a typical exponential profile 
for the atmosphere, and an aircraft altitude of 6 km, then the average refractivity, 
W, is 244 N-units for a slant range of 10 km to the ground, 249 W-units at 100 kn, 
and 255 N-units at 150 km; this error is a bias amounting to not more than 1 part 
in 10° over the ranges of interest. Therefore, eqn (9) may be calculated just once at 
some arbitrary 9, and the resulting n used for all ranges. The value of 6 that 
leads to the simplest implimentation is 90° in which case R, is equal toh = ha 


With these simplifications, and using eqn (1), eqn (9) can be integrated easily to 


obtain 
-6 6h wh - 
= 10 De mf{-6h -6h We om [{ -wh -Wh 
= = 2 + S 
n-l Gi, - %) =F, I 8 ( 8-e a) é (e 8-e a)| (10) 


It is noted that all of the parameters in eqn (10) are readily determined (or 
assumed known) except for h, and ha: These two parameters can be determined from 
the data reduction process in the system computer. After the first two or three 
passes in the data reduction computer program, the x, y, z coordinates for each 
aircraft and ground station position is known relative to a plane, the origin of 
which is located at some altitude Lo relative to sea level. The geometry is 


shown in Figure C-2. 


C-5 


ye! 
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Figure C-2. Cross-section of earth showing coordinates of aircraft and ground 
station relative to a plane, the origin of which is located at an 
altitude of Z a above sea level. 


It is seen that 


h, = d, -a 
(ae 2 2, 2 Z 
and d, = ( a7+¥)) + (2, + 2, + a) 
Likewise, (11) 
h = 
8 


d-a 
s 
x? + se + (2, tas a)? 


Eqns (11) must be calculated once and stored for each aircraft position and each 


2 
and d, 


ground station location; this requires only 50 to 70 computer words plus the 
storage of Z, and a. 

With the calculations of eqn (11) completed, the average index is found 
immediately from eqn (10), and each range correction is found using eqn (S). 
SUMMARY 


The correction procedure described above is exceedingly easy to apply. The 


following steps are. required in the computer program: 


After the ground station and aircraft positions are correctly determined 


(except for propagation errors), eqns (11) are solved for each h, and h 
and stored. 


8 


For each R(J, X) in the computer program, eqn (10) is solved for 7 - 1. 
The range correction is determined using eqn (5). 

Each range correction is subtracted from each R(/, X) to obtain the true 
range. 

Another pass in the data reduction process is initiated tn determine the 


true aircraft and ground station positions. 


The over-all accuracy of the range correction procedure has been shown to be 


approximately two parts in 10° with good radiosonde data (1 part in 10° error 


due to modeling errors, and 1 part in 10° error due to the approximation in 7). 


If the refractivity:at some arbitrary altitude must be estimated, then an error 


of 1 part in 10° for each 10% error in the estimation will also be present. 
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Appendix D 


Scientific Field Data Reduction Program 


1.0 INTRODUCTION 

The flow diagram of the position calculation program is shown 
in Figure D-1. The program may be separated into five major tasks 
namely: 

a. Parameter and data input which is processed to provide 


program control parameters and initialization data. 


local XYZ coordinates, 

c. Batch data reduction which separates the data into batches 
and locates the local XYZ coordinates of the position 
reference set and position sets along with an error analysis 
of the calculated position. 

d. Coordinate conversion to change the calculated XYZ coordinates 
into geographic lat-long or UTM. 


e. Calculate any desired auxiliary positions. 


The required input data is taken from the mission derived tape 
containing mission initialization and raw data from the position 
reference set. 

Program initialization data may be overridden by input data 


cards. Also, program flow may be modified by proper choice of 


control parameters. 


: b. Coordinate conversion to change all input position data to 
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2.0 PARAMETER AND DATA INPUT 


In order to initialize the data reduction program certain input 
data is required.’ The input tape will provide | 
a. Initial position set position estimates in either lat-long | 
or UTM coordinates plus altitude in meters. 
b. One sigma position estimates of location error in meters. 
c. Flight parameters such as turning points or actual position 
estimates and max velocity. 
d. Number of position sets. 
e. Meterological data for propagation correction. 


f. Raw data received from the position reference sets. 


Proper arrays are set up in global memory to contain this 
information so that it can be accessed by the various data processing 
routines, The raw data is unpacked and stored in arrays according 
to category such as aircraft altitude, time tags, lost lock flags, 
meterological data, auxiliary point data, message data and range 
change data. Before storing the range change data, the coarse and 
fine words are combined and normalized into one 32 bit floating 


point word. 


3.0 COORDINATE CONVERSION 

All position input information is changed into local XYZ 
coordinates. This includes the initial position set estimates as 
well as all required aircraft initialization positions that may 3 
be called out by the program. The coordinate conversion is i 


described fully in Appendix B. 
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4.0 BATCH DATA REDUCTION 


The flow diagram for the Batch Data reduction is given in Figure 


D-2 pages 1 through 8, with a detailed flow diagram of the Sub- 
routine REDUCE shown in Figure D-3 pages 1 thru 7. This is the 


major part of the data reduction task. 


4.1 FLOW DIAGRAM DESCRIPTION 

The input range change numbers represent one way propagation 
time differences from the initial measurement. These numbers are 
examined for blanks representing data drop outs and parity errors. 
This information is utilized along with the lost lock flag data to 
arrive at a lost data matrix where 1 represents a good data point 
and 0 indicates a blank data point. All bad data is set to zero 
value. The data point preceding a lost lock flag is considered 
bad in arriving at the lost data matrix. 

A test of the lost data matrix is made to determine the maximum 
number of aircraft location positions KMAX. The data is compacted 
to close up any blanks created by deletion of aircraft positions 
not satisfying the minimum station visibility criteria. An average 


propagation correction factor is utilized to arrive at an average 


propagation velocity. This velocity is used to scale the range change 


times to range change distances in meters. The time tag vector is 
converted to real time utilizing the input A/C times versus time 
tag number. The range change measurements are tested against a 
maximum number derived from multiplying the aircraft velocity 
times the measurement time difference. Any range change distances 


(RCM) exceeding this maximum are deleted and the lost data matrix 


and KMAX updated to include this deletion. 
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Figure D-2., 


Batch Data Reduction Diagram 
(Page 1 of 8) 
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Utilizing the lost data matrix, ground station locations and 


position error estimates, the data reduction is partitioned into 
batches for data reduction. The first batch is selected by choosing 
the prime stations along with the closest geographically scattered 
locations having the most good data. | ; 

The batch input data is transferred to the reduction working 
storage area and the data pass control PASS and propagation correc- 
tion indicator PSET are set to zero. PASS is incremented and the 
pass dependent program control variables are set. The integer 


program control variables are. 


Variable Control 

BOGMAX Bogus data edit iteration control 

IGC Position set location estimation 

IAC Aircraft location estimation 

IPQ Quantity of progress printout 

IAl Estimate PS oscillator frequency 

IA2 Estimate PS oscillator drift 

IPO Pass dependent print control 

NMAX Aircraft number of dimensions 
estimated 

NGAX Position set number of dimensions 
estimated 

NSTA Required minimum number of ground 


stations visable at each A/C location 


IWS Unknown station sigma reset control 
PASMAX Maximum number of data reduction passes 
DRMAX i. Was Giz? between eonmntes and actual 
RCM data bogus data editing 
ICONVR Reset convergence criteria 
D-20 
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Variable Control 

IBOGIE Bogus data editing 

ITMAX Maximum iterations of REDUCE without 
convergence 


The error estimates are set. These are as follows: 


Variable 

SIGR Range Measurement error 
SIGRO Initial range est error 

SIGA Initial freq est error 

SIGB Initial freq drift est error. 


A test is made on IAC and if aircraft positions are not being 
estimated, the aircraft estimated positions are set equal to the 
current locations. The parameter IWS controls resetting the posi- 
tion set variances. The maximum step sizes, convergence and di- 
vergence test limits are set and ground station parameters set if 
PASS = 1. The constraint weights are computed. On the first pass 
of the first batch, first try aircraft flight positions are placed 
in array Q and QE. On the second pass, the estimated positions 
QE is equated to the present computed position. 


After aircraft initialization, the following control parameters 


are set. 
ITER Iteration number 
INAL Error analysis 
IACS A/C convergence 
IGCS PS convergence 
IDIV REDUCE divergence 


t 


Chelsea iinet mem aa oh inddeia thine penta aad 


«aw sas 


The iteration counter ITER is incremented and the range change 
model is computed. The error vector DR is computed from the 
difference between the range change measurements ard computed range 
change values. The main data reduction subroutine is entered. 
Three control variables areset in this subroutine which control ; 
program flow upon exit from the routine. After first pass through | 
REDUCE, with proper convergence the error analysis IANAL parameter 

is set to 1 which causes the error analysis to be computed. If 

ICONVR is set to 1, the convergence criteria may be reset before 

calling REDUCE a second time. With successful exit from REDUCE the 

bogus data control IBOGIE is tested. If bogus data editing is 

directed, this task is performed resetting last data matrix if 

required and reiterating through REDUCE. After successful data 

eddting, a test of propagation correction control variable PSET 

and maximum pass PASMAX is made. If PASS is equal to PASMAX and - 
PSET is equal to 0, the propagation correction is made and the 

program recycled thru the last pass to refine the positions. 

Successful convergence of REDUCE with no detected bogus data 

points causes the program to return to the beginning to pick up 

the next batch for processing. For all batches after the first 


batch, the aircraft positions are held fixed by setting IAC to 0. 


4.2 SUBROUTINE REDUCE 
The flow diagram of REDUCE is shown in Figure D-3 pages 1 
through 7. For a detailed description of the process involved 


refer to Appendix A. 
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4.3 PROPAGATION CORRECTIONS 


Appendix C details the propagation correction methods used 


in the program. 


4.4 COORDINATE CONVERSION 
The coordinate conversion routine is reentered to convert the 


local XYZ coordinates to Lat-Long and UTM before printing out the 


results of each batch. 


5.0 AUXILIARY POSITION CALCULATION 

After successful completion of the data reduction a test is made 
to determine if any aux positions are desired. If some positions 
are required, the position set location requesting the position is 
used as the origin of a local XYZ coordinate system. In other 
words the point X = 0, Y = - and Z = O correspond to the Lat-Long 
and Elevation of the PS location. The X Y Z of the AUX point is 
computed using the range, with elevation and azimuth angles as shown 


in Figure D-4. 
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Figure D-4. Aux Point Local Coordinate Calculation 


These local coordinates are then converted to Lat-Long or UTM 


by the coordinate conversion routine. 
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APPENDIX E 


FIELD DATA REDUCTION PROGRAM 


With the exception of operator controlled program flow and 
bogus data editing, the field data reduction program is the same 
as the scientific data reduction program described in Appendix D. 

A detailed flow diagram of the field data reduction portion 
reflecting these modifications is shown in Figure E-1l, pages 1 


through 8, 
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APPENDIX F-1 . 4 
I LRPDS RELIABILITY MATHEMATICAL MODEL DEVELOPMENT 
1. Assumptions 
The following assumptions were applied in the development of the 
mathewatical model. 
a. The system is a series system as described by the LRPDS Reliability 
Block Diagram (See section 10-2 of this report.) y 
b. Operating times and duty cycles for each subsystem and unit used 
during a mission are as presented in Table F-1.2._ 
c. Times to failure follow an exponential distribution. 
d. "Reliability" is defined es the probability that an item will 
operate without failure for the duration of the mission. 
2. LRPDS Mission Time 
The LRPNS consists of three subsystems, each of which operates for 
a different period of time during the mission. Furthermore, each 
subsystem is composed of units operating for different periods of 
time, Therefore, individual "mission" time has been developed for 
each subdivision of the system. These time elements are based on the 
“standard mission" as defined in paragraph 10.2.1 of this report, and 
are presented in the Reliability Block Diagram. 
The operating times for each system subdivision were derived based 
on mission timing data devdloped during the LRPDS Mission Description 
study. The results of this study hav? boen usod to dotermino tho 
approximate duration of cach element of a “standard mission" (See ; 


table F-1.1.) These time elements are summed as appropriate to 


establish mission times for each unit of the system (See table F-1,2). 


TABLE F-1.1 V 
TIME ELEMENTS OF STANDARD MISSION 


se MISSION ELEMENT TIME (Minutes) 


Oscillator Warm-Up 


PS's and PCC/PS 240 
A/C 120 


Mission Initialization, etc. 


PCC 30 

A/C 30 

PS's 5 
Initial A/C Position Run 7 
Standby for P2 Command 5 
Ranging Mission 46 
Message Exchange 9 
Data Processing & Reduction 15 2s 


Estimate of computer time required after completion of 
mission, 


TABLE F-1.2 


I — 


EQUIPMENT OPERATING TIME 


* OPERATING TIME 
SUBSYSTEM/EQUI PMENT DUTY CYCLE (HOURS) 


7. PCC Subsystem 


* Master Oscillator 1,00 5.9 

j PCC/PS (All except XMTR) 1.00 1.6 

¥ PCC/PS Transmitter 03 1.6 
All other PCC Units 1,00 1.9 | 


Reference Position Subsystem (A/C) 


Rubidium Frequency Standard 1.00 3.5 
Transmitter 75 1.6 
All other A/C Units 1,00 1.6 


Positioning Set Subsystem 
Master Oscillators 1.00 5.2 
I Transmitters 02 1.2 
All other PS Units 1.00 1.2 


In addition, operating duty cycles were determined for those items 
(transmitter) that are operated intermittently. These are based on 
the approximate time required to transmit each of the various commands 


and messages associated with the standard mission. 


3. Basic Reliability Model 
The reliability of a system composed of a number of elements 
connected in series (i.e., a failure in anyone of the elements causes 


a system failure) is given by the expression: 
N 


Rect) * JA a Pee) Q) 


Where: 


Ro(T) = the reliability of the system for mission time T. 


Ri (ti) = the reliability of system element i for its mission. 


time ty tyST 
N = the number of elements in the system 
This basic model can be expanded by considering an individual 
element of the system. Assuming the exponential distribution, 


Rt) = exp bak t) ; (2) 


Where: 
X= the failure rate of the system element under the opcrational 


and environmental conditions encountered during time period 


t. 


a eal = SS SSS 


ae 


— 


— 


Expression (2) assumes that the environmental and operating conditions 
remain constant for the entire mission. In gencral, individual 

subsystems of the LRPDS can be assumed to operate in & constant 

environment for the duration of the mission, seinen, certain system 

elements, in particular transmitters, will be cycled on and off during 

the mission. Assuming the failure rate of an item during operation 


is different than during periods of non-operation in the same environ- | 


ment, expression (2) can be modified: 


Ree) = exp {-LXolt-(d/e)) #Xg_g (tC = ae) 3) 


n-o 


Where: 
d/c = the duty cycle of the time, where duty cycle is the ratio 
of operating time to total time. 
Xo = the operating failure rate 


>.n-o = the non-operating failure rate 


Failure rates arc neown diy determined under operating conditions 
and non-operating failure rates are not available in most cases. 
Therefore, expression (3) is not readily useable. However, a number 
of non-operating or storage life studies have indicated that operating 
failure rates are at least 10 times greater than non-operating failure 


rates in the same environment. Therofore, for the purposes of this 


1. 


Onc such study is reported in RADC-TR-67-307, “Dormant Oporating j 
and Storage Fffects on Electronic Equipment and Part Roliability". 


model, it will be assumed that A gt oh Age in which case, expression 


(3) becomes: 


Ree) 7 exP fl At(d/c) + .1 At CL - a/e)1} 


Where: 


d= the operating failure rate. 


Rearranging this expression gives: 


Rc) 7 exp [-At (.9 (a/c) + .1)) (4) 


Expression (4) is used as the basic elemental expression in the 


development of the LRPDS Reliability Model. 


4. Detailed Model Development 
The LRPDS Reliability Model is based on the reliability block A 

diagram (Figures 10-1 through 10-7 of this report). The system is_ as 

assumed to have the following configuration: 

a. One Position Computing Central (Block number 1000 of the block 
diagram). 

b. One Reference Position Set (Block #2000) 

c. 27 Positioning Sets (Block #3000), only 5 of which include a ‘ 
Data Display Unit (Block #3300). 


From expression (1), the system reliability is expressed as: 


Re(5.9) ~ ®1000(5.9) * ®2000(3.5) ° ®3000¢5.2) (5) y 


71 
The 4-digit numbers refer to block numbers of the reliability 


a rr 


block diagram, and the numbers in parentheses are the mission times, 
in hours, for the respective subsystem. “° 

The expression R3000(5.2) referes to the entire PS subsystem 
consisting of 27 Positioning Sets. Each of the subsystems reliabilities 
is expressed as the product of the reliabilities of the individual = 
units of the subsystem. Individual factors of these expressions are 
then factored into more elemental ‘subdivisions until the subsystem 
reliability can be expressed as the product of a set of exponential 


expressions each having the form of expression(4} This development 


zor the three subsystems is described in the following paragraphs. 


4.1 Position Computing Central Reliability Model: 


Rioo0 ~ F1100 - Fiz200 - i300: ®1400 


(Note, block 1500 represents GFE and is not being 
considered in this interim report. Therefore, Ri 800 and 


Ri900 is set equal to 1 in this model.) 


Ri 920 can be expanded according to Figure 10-3, giving: 


Ry220 = 21201 - Ria22 - Fize3 - Piso 


The PCC model has now been factored into the lowest terms possible 


at the present state of system definition. Each factor is now expressed 


2. This subscribing method will be continued without further explana- 


tion in the remainder of this appendix. 


rg Oe eT OP 


in the 


values 


R991" 
Ri 2227 
Ri 223 


Ri 220 


R1300 
Ri 400 
R500 


Ri 230 


exponential form by substituting appropriate time and duty cycle 


from Figures 10-2 and 10-3 in expression (4): 


exp [- pe (1.6 x 1)]) = exp (1.6 Aj 591) 


en {-Aj 209. [1.6 (.9 x .03 + -1))} = exp (-.2A, 990) 


= exp (-1.6 digas) 


exp (-5.9 A, 990) 


exp (-1.9 Pade) 


exp (-1.9 Ly 


exp (-5.9 A 1336) 


Substituting these values in the expression for Ri 000° collecting 


terms and simplying: 


4.2 Reference Position Set Reliability Model; 


Roooo * F2100 ° 2200 ° F2300 ° ®2400 


R9200 


R 


+ 5.9 (A 1920 + 


1000 = &xP [- [1-6 CAjo97 + Ayea3) * +2 Azooe 


) + 1.9 CA 


: \ 
1224 * A1500 1300 * 41400) || 


4 


can be expanded according to Figure 10-5. 


Roo00 ~ Ro210 * Re220 ° Reese + Re240 


a Roo10 ~ Fe211 + Rea12 


Expressing each factor in the exponential form by substituting 


appropriate time and duty cycle values from Figures 10-4 and 10-5 in 


expression (4): 
Roi00 7 eXP (-1-6 Ao190) 


Rooi1 7 ©xP (-1.6 Ago])) 


- \ - ie ” Ss i} 
Roo12 exp | ho012 [1.6 (.9 x .75 + -1)]} exp (-1.2 9912) tj 


Roga0 = eXP (-1-6 A g990) 


R2230 ~ 2230) ete q 


| Rooag 7 XP (1-6 Agogo) 
Roosq * xP (-1.6 Agoso) 


F 
E 


Roggg = XP (-1.6 Agggo) 


Substituting these values in the expression for Ro go00° collecting 
terms and simplifying: 


Rogoo 7 xP {7 [2-6 €Aoi09 *A2211 *A 2290 * 2230 * Az240 * Azeso * 


A + 1.22 


2400 2212 * 3-5 Agzo9 || 


4.3 Positioning Set Reliability Model 
The Positioning Set subsystem includes 27 PS's, only 5 of which 
ts include DDU's. Therefore, the subsystem reliability model is expressed: 


; 27 5 
R3o00 * ‘83100 ° ®s200 * 83400? ° ‘3300? 


as da ats emai e a  et lla  ae ea e  ataln 


ee 


nnn 


R 


can be expanded according to Figure 10-7. 


3100 


R3100 ~ 83110 - ®3120 ° ®3130 


R3110 ~ 83111 ° ¥3112 


Expressing each factor in the exponential form by substituting | 
appropriate time and duty cycle values from Figures 10-6 and 10-7. 
in expression (4): 


Ray71 7 exp (-1.22377)) 


R312 7 exP |- Aar12 [1.2 (.9 x .02 + -1) ]| = exp (-.14A5))9) 


R320 = exP (-1-2A 3190) 
R3130 ~ exP (-1.2A3) 30) 
R3200 = exP (-5.2 3900) 


R (-1.2A 


3300 ~ ©&*P 3300? 


R3400 = 1 tor this model because battery is not yet defined. 


Substituting these values in the expression for R30900° collecting 


terms, and simplifying: 


R3000 ~ {exp[-1.2 (A3111 + 3120 +4 3130) ~ +142 sirz 


27 
3200 | | 


: (enpi=1.2 3300) | 


-5.2 A 


5. Relationship Between Elements of Reliability Model and Hardware 


Subdivisions 


The failure rate factors ty represent failure rates of the 


xxxx? 
largest subdivisions of major units of the respective equipment for 


F-10 


\ 


ee | ,=— ream UlUlC Ell eles 
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which the mission time and duty cycle would be consistent, and corres- 
pond to the various blocks of the reliability block diagram. In 
essentially all cases, these blocks can be further subdivided into 
modules; however, such subdivision is not necessary in most cases 
because the block failure rate is equal to the sum of the failure rates 
of the modules. However, in the case of the RTU's, modules making 
up the transmit circuitry operate at a lower duty cycle than modules 
making up the receive circuitry. Module or circuitry operating at 
the different duty cycles are listed below. 
Receive Circuitry (Blocksl111, 2211 and 3111) 

Power Converter 

RF/IF 

Synthesizer (Receive portion) 

Carrier IQ/VCO 

Code Detector 

CACU (70% of total module failure rate) 

Clock IQ/Phase Control 7 

RTU Cable Harness 
Transmit Circuitry (Blocks 1112, 2212 and 3112) 

Power Amplifier/Modulator 

Synthesizer (Transmit Portion) 


CACU (30% of total module failure rate) 
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APPENDIX F-2 ss 
LRPDS RELIABILITY ALLOCATION 


1. Allocation to Subsystems At 
The effective maximum allowable system failure rate is 26000 


1. This failure rate is 


failures per million hours (26,000 £/10°), 
allocated to the three subsystems such that the sum of the failure 
rates of the individual subsystems is equal to the. total system 
failure rate, Other factors considered in the allocation include the 
relative subsystem complexity and variation in mission time constraints. 
Two basic relationships are used in performing the allocations. 
First, in a series system, the system reliability is the product of 
the subsystem reliabilities. Assuming an exponential distribution, 
N 
exp (-A _T) = Ul exp (-A,$,) exp [- Yt] : 
i=l i-1 ' 
Where: 
Me = the ayet ee failure rate 
T = the system mission time (the maximum mission time 
for any subsystem) 
Aa = the failure rate for sybsystem i 
ti = the mission time for subsystem i 
N = the number of subsystems 


Taking the logarithm of this expression and rearranging gives: 


N 
r~ = Asts (1) 


x The failure rate of a series system is the reciprocal of the MTBF. 
Thus, the effective maximum allowable failure rate for the LRPDS = 
1/38 hours = ,026 failures per hour. 


F-12 [| 
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Furthermore, in a series system the system failure rate is the sum of 


the subsystem failure rates such that 


N 
- ‘ (2) 
Xs se As 


Both expressions (1) and (2) are used in performing the allocation. 
The system and subsystem mission times being assumed for this 


report arc as follows (sec paragraph 2.1.2). 


LRPDS System: T = 5.9 hours 

PCC Subsystem: ti 000 = 5.9 hours 
RPS Subsystem: toooo = 3.5 hours 
PS Subsystem: t3000 = 5.2 hours 


Substituting these values, and the system 


(A, = ,026 f/hour) 


failure rate (AQ = .025 f/hour) in expressions (1) and (2): 


5-9. j000 , GeS aon i S28 Asean 
5.9 5.9 5.9 


= .026 f/hour 


' ) - 


Ax000 + 42000 + A3000 = -926 (4) 
These two expressions cannot be solved simultaneously because (1) 
there are 3 unknowns, and (2) direct solution wouid provide at least 
one negative A value. Therefore, x valucs are determinod that will 


satisfy expression (3) and that reflect tho relative subsystom 
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complexity. Different environmental "K-factors" are also considered. 
Final allocations are determined by adjusting the nN values to satisfy 
expression (4) while maintaining the relative weights that were 
established using cxpression (3). 

The following observations are made regarding the complexity and 
other factors relating to the relative failure rates of the three 
subsystems. 

a. Comparison between RPS and PS: The RPS isnot more than 10% 

more complex than individual PS. Therefore, an individual PS 

failure rate (,\ 3000 would be approximately 3. sacs in the 

same environment. However, the PS subsystem consists of 27 

individual PS's, Thus, the PS subsystem is estimated to be 

27 x .9 = 24,3 times as complex as the RPS subsystem. Also, 

even though the two subsystems operate in different environments, 
the environmental "K-factors" are the same for both subsystems, 2° 
Based on these observations, the RPS and PS subsystem failure 


rates are allocated such that: 


\ a ) 
AS 3000 24.3," 9000 (5) 


b. Comparison between RPS and PCC: Based on the number and types 
of blocks in the PCC block diagram (Figure 10-2) as compared to the 
RPS block diagram (Figure 10-4), it is apparent that the PCC is 


at least 5 times as complex as the RPS. However, the PCC operates 


See section 10.4 for a discussion of reliability "K=-factors" 
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in a fixed, ground environment while the RPS operates in an airborne 
environment. The average ratio between fixed ground and airborne 
"K-factors" in MIL-HDBK-217A is approximately .16:1. Based on these 


observations, the system failure rates will be allocated such that: 


Arooo * (5 * -16)A so90 


rw 


Manos “ 2000 ; (6) 


Substituting (5) and (6) in (3): 


} a 
dgaag + 50. was + 24.3(.88) oo99 = 026 


2926 . (9011 £/hour (7) 


A = 
2000 ~ 0 9 


Substituting (7) in (5) and (6): 


A 3000 7 24-3 (0011) = .0267 f/hour (8) 
A,o00 = °8 (+0011) = .0009 £/hour (9) 
These values are adjusted such that expression (4)' is 
satisfied: 
A 1000 * 2000 +3000 = +0009 + 0011 + .0267 = .0287 


But, from (4): 


A,000 +2000 * A s000 7 +026 


.026 
» 0287 
slightly to provide the following allocations. 


Therefore, the >. values are multipliod by = .91 and adjuated 
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e tJ = 6 
Pec: Ajoo9 = +0008 £/nour = 800 £/10 


rps: 2000 7 °0009 £/nour = 900 £/10° 


6 
ps: 3000 7 °0243 f/hour = 24,300 £/10 


Total System \ = 26,000 £/10° 


‘2, Positioning Set Subsystem Allocation : x 

The PS subsystem failure rate (Nena) is allocated to the i 
lower subdivision first because the Positioning Unit is a fundamental 
element of the system, and similar units are used in the RPS and PCC 
subsystems, 

The PS subsystem consists of 27 Positioning Sets, 5 of which 
include a Data Display Unit (DDU). Thus, the subsystem failure rate 
is allocated to 27 positioning units (block 3100 of Figure 10-6), 
27 Crystal Oscillators (block 3200) and 5 DDU's (block 3300). The 
battery (block 3400) is ignored in the initial allocation because (a) ue 
specific type of battery to be used is in doubt at this time and 
(b) batteries are highly reliable devices (typical failure rates are 
in the order of 1 £/10° or less) and therefore, do not contribute— — 
significantly to the system failure rate, In the final allocation, a 
failure rate of 1 £/10° is allocated to the battery. 

The PS subsystem reliability is expressed in the model (paragraph 


10.2.3c) as: 


clio .48t 5 
R3300 * “83100 ° ®3200? ‘33007 (10) 
Where: 
R3100 7 xP | - (1.2 A399) (11) 
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R3000 = cxp [ -(5,2 fren (12) 
Rg309 ~ xP [- hie Agana)! (13) 


Also, the reliability block diagram (Figure 10-1) gives Tojo © 
5.2 hours, and the allocated failure rate for this subsystem is 


A aie = ,0243 f£/hour. 


Therefore: 


R3y99 ~ exP [- (.0243 x 5.2)] (14) 


Substituting (11), (12), (13) and (14) in (10): 


exp [ -(.0243 x 5,2] = exp [ -(27 x 1.2) Asig9 + 


Taking the logarithm and rearranging: 


) } we 
19.8 A aus #27 A opgg t 1218 Aeggg = -0249 (15) 


Also, this is a series system and, 


) = 
27 d3199 + 27 y\ soot 8 d ee66 0243 (16) 


‘The following observations are made concerning relative com- 
Pplexitics. 

a. Comparison of Crystal Oscillator and Positioning Unit, The crystal 
oscillator consists of approximately 150 parts compared to 
approximatcly 2500 parts in the Positioning Unit. Therefore, it 

is assumed that: 


3000 = sean’ = .06A (17) 


3100 
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b. Comparison of DDU and Positioning Unit: The DDU contains 
r approximately 140 parts, 20 of which are LED indicators. These 
devices are roughly equal to 14 diodes and one IC, or the 
equivalent of 16 parts. Therefore, it is assumed that the DDU is 
approximately 500/2500 as complex as the Positioning Unit, and: 
43300 = +2A 3100 (18) 
Substituting (17) and (18) in (15): . 
19.81.5199 + 27 (.06Agi 99) + 1618 (.2 Agipo) = -0243 
-0243  _ 
3100 ~ Ses. 7 100223 | 
Substituting this value in (17) and (18): 
t Agze09 7 +08 x -00113 = .00007 
= = \ 
rg309 7 2 * -00113 = .000226 a 
This is checked in (15) which states: 
27) + 27) 475A = .0243 . 
“3100 3200 3300 : 
But, substituting the calculated ae values; & 
27 x .00113 + 27 x .00007 + 5 x .000225 = .0335 
2 -0243 3 
Therefore, the ’. values are adjusted by a factor of “5555 Ca ach 
giving: 
: 6 ‘ 
Positioning Unit: ). 3199 7820 £/10 
nN 
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Crystal Oscillator: anne = 50 2/10° 


6 
DDU: ? rs = 169 £/10 


Battery: Aaaee wl £/10° F 
2.1 Positioning Unit Allocation 


The positioning unit failure rate, .\ 5,99 ~ 820 £/10° is allocated 
to the RTU, DPU and control panel in proportion to the eStimated 
failure rates as determined during preliminary predictions. This 
procedure is used in lieu of a complexity or part count weighting 
procedure because the analog circuitry of the RTU is expected to 
have considerable higher failure rate to part count ratio than the 
digital circuitry of the DPU. 

A preliminary prediction that was available before this allocation 


was performed indicating approximate failure rates as follows:° 


Unit Failure Rate Percentage 
RTU 680 ~60 
DPU 210 : 20 
Cont. Panel 220 ; 20 
Total 1110 1.00 


3. These are preliminary and are used for weighting quantities only, 
See section 10.4 for actual prediction results. : 
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Using these percentage as weighting factors provides the following 


allocation of Positioning unit failure rates: 


RTU: = 820 x .60 = 490 £/10° 


A 3110 


DPU: = 820 x .20 = 165 £/10° 


r 3120 


Gant. Panel. > = 820 x .20 = 165 £/10° 


3130 


2.2.1 Allocation to Modules of PS Positioning Unit: 
The RTU and DPU allocations are allocated to modules using 
approximate failure rates obtained from the preliminary predictions. 


These allocations are listed in Table F=-2.1. 


3. Reference Position Set Allocation 


The Reference Position Set consists of a Control and Monitor Unit 


“(block 2100 of Figure 10-4), a Reference Position Unit (block 2200), a 


Rubidium Frequency Standard (block 2300) and an Altimeter (block 2400). 
The subsystem reliability can be expressed according to the model (para- 


graph 2.1.5b) as: 


= Rh R 


Ro000 ~ 82100 * R2200 + ®2e300 ° F2400 (19) 
Where: 

Roi00 7 exP (-1.6A 5199) (20) 

Roo90 7 xP (-1-6A g999) (21) 

Ro3909 = xP (-3.5A 4590) (22) 


Rogoo = &xP (-1-5 Aggoo) 
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Also, the reliability block diagram (Figure 10-1) gives Too00 = 3.5 
hours, and the allocated failure rate for this subsystem is 


6 
r s000 7 900 £/10°= .0009 £/hour 


Therefore: 


R = exp [- (.0009 x 3.5)] (23) 


2000 
Substituting (20), (21), (22),and (23) in (19): 


exp [= (.0009 x 3.5)} = exp [- (1.6A 3199 + LCA sang + 


) 


Taking the logarithm and rearranging: 


) . 
-46 (9100 +/A2200 * A2400) * 2300 = -9°°9 (24) 

Also, the RPS is a series system, and 
A100 +4 e200 +4 2300 * A200 = +9008 (25) 


The following observations are made regarding therelative 
complexity of RPS units. 

a. The Reference Position Unit is similar to, but somewhat more complex 
than the PS Positioning Unit. The unit contains more digital 
circuitry, a more complex memory unit, and a separate memory 
power supply. However, the unit does not contain the 33 data 
switches such as are included on the PS control panel. Therefore, 
the Reference Position Unit should be only slightly less reliable 
than the PS Positioning Unit. The preliminary allocation for the 


positioning equipment is; 
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r / ne = ORES f/hour .. (26) 
Nw 
b. The CMU is a relatively simple control box containing switches 
and indicators for operating the RPS. The unit contains 
approximately 7 switches, 11 Bite indicators, and 2 electrical 
connectors. Failure rate values for this type of component is 
usually in the order of 1 to 5 £/10°, so that the total CMU 
failure rate should be between 20 and 100 £/10°, Therefore, the 
preliminary value for the CMU failure rate allocation is: 
Ao100 * .0005 f/hour (27) 
c. The altimeter is a relatively simple device consisting of a 
transducer and simple A-D converter. The failure rate of this 
device should be significantly less than that of the CMU. Therefore, 
the preliminary value assigned to this unit is: 
se Ao400 7% - 900005 f/hour (28) 
Substituting (26), (27) and (28) in (24): 
.46 (.00150 + .00005 + .000005) + Ag300 = .0009 
A9300 * .0009 - .0007 = .00025 f/hour (29) 
The sum of the values of (26), (27), (28) and (29) is 
,00150 + .00005 + .000005 + .0002 = .001755 
But, from (25), this sum should be .0009, 
Therefore, the preliminary A} values are adjusted by a factor of 
uO? a .00, and the allocated values are: 
.0018 
6 
CMU: A 2100 ™ 23 £/10 
6 
Ref Pos Unit: Agggg * 750 £/10 
6 
wK Rubidum std: Aggqq * 125 £/10 


Altimeter: d 2400 » 2 £/108 
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3.1 Reference Position Unit: 
The Reference Position Unit (block 2200) consists of five units 
as shown in Figure 2-5. The RTU is identical to the PS RTU and is 


allocated the same failure rate, Such that: 


Maate = ,00049 f/hour': 


The DPU (block 2220) is similar to the PS DPU in construction. 
The memory, which is included as part of the PS DPU, is a separate. 
unit in the RPS. However, there are approximately twice as many 
digital modules, each of which is similar in complexity to the PS 
DPU modules. Therefore, the preliminary allocation for the RPS DPU 


is approximately the same as for the PS DPU, or: 


epee = .00015 f£/hour 


The remaining failure rate of .00019 is allocated to the memory 


(block 2230), memory power supply (block 2240), and the Antenna 
Filter (block 2250) in relation to the estimated relative complexity. 


The preliminary allocation for the positioning equipment is: 


RTU: p - be Gi 490. £/10° 
DPU: a d e999 7 150 £/10° 
Memory: A wang m 75 £/10° 
Memory P.S.: ” gag = 33 2/109 


: ‘ ) ‘ 6 , 
Antenna Filter: “9959 2 £/10 
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3.1.1 Allocation to Modules of the Reference Position Unit 

Failure rates of the RTU are allocated to modules in accordance 
with the allocations to comparable modules of the PS RTU (block 3110), 
and are summarizcd in Table Bel. (Allocations to low subdivisions 


of the othcr units will be performed at a later date.) 


4. Position Computing Central Allocation 

Due to the preliminary status of the PCC subsystem definition, 
the allocation will be performed for the PCC/PS (block 1100) only. 
The remaining portion of the subsystem failure rate will be allocated 
to the other units when the subsystem is sufficiently defined. 

The PCC/PS is essentially the same as the PS, (including position- 

ing unit and master oscillator). However, the sheltered ground environ- 
ment, and the environmental "K=factor" is less by a factor of approxi~ 


mately .16:1. Also, the PS control panel is not used. 


Therefore, the allocation to the PCC/PS is .16 times the PS, or: 


Ar220 = .as (Asieo +A 000A = .16 (820 + 50 - 165) 


rAx229 ~ 115 2/10° 


This is allocated to the various units and modules of the PCC/PS 
by multiplying the respective PS module failure rates by .16. The 


results of this allocation is presented in Table F=2,1, 
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| TABLE Fr-2.\ FAILURE RATE ALLOCATION SUMMARY 
| SYST ELA Block | FAILURE RATE ALLOCATIONS (f/)U4HRS) 
7 REALL OA NUHBER SYSTEM SUFSHS MAIR 
I AN/US@-5¢(V) Tn a ee aa 
AN/TS9 -100 1000 ; 
Pcc/Ps 1220 ne 
| RTV . ¥o ; 
POWER CONV. — 3 
RF/IF - + 
i PA/MOOULATOR - 1s 
FREQ. SYNTH(TEAR —- Le 
CARR 1a/ VCO - 9 | 
i CODE DETECTOR - 3 
CACU - 4 
i CLOCK 1Q/PH.CONT | — = 
DPV 1223 27 
“Re BY | -- 4 
[ POWER SUPPLY | = 7 
WoRD CONTROL - a 
MESSAGE OUTPUT | = ae 
| DATA ASS ZITGLER - % 
COMMAND DECODER — 2 
MEMORY - 3 
| CABLE HARKNESS cad 4 
CRYSTAL FREQ. SIP. * [1224 
f ALL OTHER PCC UNITS b &s~ 
! 
AN/ ASO-148 | 2000 
| CONT.- MON|TOR UNIT ; 2/100 23 
REFERENCE POSITION UNIT 2200 750 
RYv red 
fi POWER CONVERT. | = 
RF/I\F - 
| PA/MAPVYLATOR = 
FRE@.SYNIUACTES) = 
CARR IQ/VCO - 
| CODE DETECTOR | - 
CACY | - 
CLOCK IQ/HH.CONT! = 
| CABLE HARNESS - 
DPU £798 
MEMORY 2220 
| MEMOLY PS, 22490 
ANTENNA FILTER 26.350 
BUBIDIVUM ©REQ STD 2309 15 
7 . ALTIMETER 2400 a 


eens cere © ome OM Of MRR sane eee ae a oe Ome ee tee eee wo we ee 


| yeas 


SYSTEM BLOCK FAILURE RATE ALLOCATION (E/ 108 WAS 
BREAKDOWN NUMBER SYSTEM Supsvs iz a  [eonuee | 


Wa SSS a rs SR EELS LAS LTS Sh 


AN/PSQ- ‘lou 3000 
FULL SUESYS TEM (27 Poss ' 249,300 
AND 5Dpy's) 
POSITIONING UNIT. 3/00 620 
RTv °- 3110 
POWER CONV - 
RF/|F - 
PA/MOOULATOR - 
FREQ@.SYNTH(TER)| 
CARR 10/VCO - 
COCE DETECT IR = 
CACU - 
CLOCK 1Q/PH.CONT{ — 
CABLE HARNESS ~ 
DPU 3120 


RCbu = 
POWER SUFPLY ~ 
WORD CONTAOL _ 


ASG OUTPUT = 
DATA ASSENGLER = 


CMD DECODER ~ 
MEMORY ~ 
CABLE HARNESS _ 
CONTROL PANEL 3130 
CRYSTAL FRE STD 3200 50 | 
DATA DISPLAY UNIT 33085 10? 
we BATTERY 3400 ( 
{ 
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APPENDIX F=3 4 
RELIABILITY CALCULATIONS 


_ =e 2 


Aa RELIABILITY MODEL SOLUTIONS 
"Standard Mission" reliability values are calculated below for 
allocated failure reates, and for worst case, typical operation, 


and reliability test prediction failure rates. 


teal PCC MISSION RELIABILITY 
The PCC Reliability for the standard mission is calculated 
using the allocated failure rates (See Table F-2.1). 


The reliability mathematical model for the PCC subsystem is: 


C ‘ 


+di200 * Ais00 * 2-9 (Aas00 *41400? |} 


( 


A values are determined as follows from allocated failure rates, 
and considering the hardware items included in each block of the 


reliability block diagram. 


Block Item } Failure Rate (F/10°) 

1221 RTU Receive Circuitry 
Power Converter 8 
RF/IF 4 
Receive Synth ll 
Carr, IQ/VCO 9 
Code Detector 3 
70% of CACU 3 
Clock 1Q/PH. Cont een 

Total 4444 43 
1223 DPU 27 
F-27 
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Block Item Failure Rate (£/10°) 


1222 RTU Transmit Circuitry 
PA/Modulator 25 
Trans. Synth ll 
30% of CACU 1 
Total Desai 37 
1224 Crystal Freq. Std. - 378 
All other PCC Units 685 


Substituting these values in the model: 


Rioo0 = exP ft 1.6 (,00007) + .2 (.000037) + 5.9 (.000008) 


+ 1.9 (.000685) ] 


= exp (-.00147) = .9985 


1,2 RPS MISSION RELIABILITY 

The RPS reliability for the standard mission is calculated 
using allocated failure rates from Table F-2,1, and Worst case, 
typical operation, and reliability test failure rates from the 


prediction data sheets in Appendix F-4, (See Table F=3.1). 


The reliability mathematical model for the RPS subsystem is: 
= exp f 
Raooo ~ xP; -[1-6 (A 2100 tA ga12 tA 2220 +As230 + \o240 
roa5o) + 1+2Agor0 + 3.5A gg00)| 


Substituting these data in the reliability model gives allocated 


and predicted mission reliability as follows: 


F-28 


f 
y 
t 


Eee Ee FFE ST wa FB 


a. Allocated Reliability: 
R = exp (4 1.6 (.000025 + .000291 + .000150 + .000110 
2000 b 
+ .00006GO + .000020) + 1.2 (.000219) + 3.5 (000025) }} 


R = exp (.0014) = .9986 


2000 


b. Worst Case Prediction 
Ro 000 = exp Co 1.6 (.00006G1 + .000545 + .000295 + .000250 
rt 
+ .000092 + .000001) + 1.2 (.000402) + 3.5 (.001321)} 


c. Typical Operation Prediction 
f. 

R = exp }-{ 1.6 (.000024 + .000275 + .000115 + .000098 
2000 j 


5 
+ .000056 + .0000002) + 1.2 (.000202) + 3.5 (000486) }, 


d. Reliability Test Condition Prediction 


c 
Roo00 = --{ 1.6 (.000006 + .000124 + .000064 + .000092 


| 


+ .000014 + .0000001) + 1.2 (.000075) + 3.5 (000085) }; 
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TABLE F-3.1 J 
RPS FAILURE RATE VALUES * 
é 
~ A (£/10°) 
BLOCK ITEM ALLOCATION WORST CASE TYPICAL R TEST — 
2100 Control & Monitor 25 60.75 24.14 5.72 * | 
Unit 
2211 RTU Receive Circuitry : 
Power Converter 44 79.51 41.37 9.82 ! 
RF/IF 19 36,01 — 18.19 7.63 y 
Receive Synth 66 123,09 65. 98 23.38 — 
Carrier IQ/VCO 54 98.61 52.46 25.78 
Code Detector 16 30.55 13,09 13.67 
70% of CACU 14 29.24 10.18 7.28 
Clock IQ/ph cont. 66 125.43 63.50 31.80 - 
Cable Harness 12 22.26 10,22 4.35 
Total Aoo11 291 544,69 274,99 123.71 . 
2220 DPU 150 295,44 114,99 63.98 
2230 © Memory 110 250.43 98 . 37 91.58 
2240 Memory Power Supply 60 92.11 55. 68 14.18 
2250 Antenna Filter 20 1,24 18 ell 
2212 RTU Transmit Circuitry 
PA/Modulator 147 267.90 132.86 48.62 
Transmitter Synth 66 121. 58 64.47 23.06 - 
30% of CACU . 6 12.53 4.37 3,12 
Total A o919 219 402.01 201.70 74.80 
2300 Rubidium Freq. Standard 25 1321.15 486.23 85.15 - 
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3120 
3130 


3112 


3200 


3300 


ITEM A 


RTU Receive Circuitry 
Power Converter 
RF/IF 
Receiver Synth 
Carrier IQ/VCO 
Code Detector 
70% of CACU 
Clock IQ/Ph. cont 
Cable Harness 


Total fae 
DPU 
Control Panel 


RTU Transmit Circuitry 
PA/Modulator 
Transmitter Synth 
30% of CACU 


Crystal Oscillator 


Data Display Unit 


TABLE F-3,2 
PS FAILURE RATE VALUES 


LLOCATION 


50 


170 
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WORST CASE 


TYPICAL 


R TEST 
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C 1.3 PS SUBSYSTEM MISSION RELIABILITY 


as The PS subsystem reliability is calculated using allocated failure | 
rates from Table F-2.,1 and worst case, typical operation, and 
reliability test failure rate data from Appendix F-4. (See Table 


F-3.2). The reliability mathematical model for the PS subsystem is; 


27 5 
R3000 7 83100 * ®3200) ‘833002 


Gq 
R3100 * xP. -[1.2 OA 311 +A3120 +A3130) + DNgr a) 
R - exp; -[ 5.2 (A yy 
3200 Bey mE 32007?) 
( a. ) 
R3300 ~ exP, 11-2 (A 3390) }4 


The values obtained from Table F-2.1 and Appendix F-4 are listed in 


we Table F-3,2. Substituting these data in the reliability model gives 
allocated and predicted mission reliability as follows: 
a. Allocated Reliability: 
R3100 = oxp { -L 1.2 (.000279 + .000165 + .000165) + .14 (.000211) J. 
R3100 = exp (-.00076) = ,99924 
Ryoo9 7 exP | -[5.2 (.000050)}', = exp (-.00026) = .99974 
R3309 = exP tae (.00017)] : = exp (-.000204) = .9998 
Rao09 7 [¢-9992) ¢ .9997)}77 (.9998)° 
/ Raggy ~ (-9989)77 (,9998)° = .9696 


“€ 
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b. Worst Case Prediction 


R3100 


R3100 
R3200 
R3300 
R3000 


c. Typical 


R3100 
R3100 
R3200 
R3300 


R3000 


d. Reliability Test Condition Prediction 


R3300 


R3100 
R200 
R43300 


R3000 


Operation Prediction 


exp{ -[ 1.2 (.000538 + .000339 + .000257) + .14 (.000401)}) 
exp (~.00142) = .9986 
) 
exp | ~[5.2 (.000374)], = exp (-.00194) = .9981 
f » : 
exp 7 =f 1.2 (.000165)]}< = exp (-.000198) = .9998 


[(.9986) (.9981)}?" (.9998)° = .9126 


Z \ 
exp; -(1.2 (.000265 + .000132 + .000108) + .14 (000202) ]} 
exp (-.00063) = .9994 

C ay 
exp : -(5.2 (.000074)]. = exp (-.00039) = .9996 
exp C[1.2 (.000096) } ‘ - exp (.000115) = .9999 


[(.9994) (.9996)]7" (.9999)° = .9722 


>) 
exp fat 1.2( 000124 + ,000075 + .000041) + .14 (,000074)); 


exp (~.000308) = .9997 
‘y ° 
exp | -[5.2 (,000027)]: = exp (-.00016) = .9998 


exp T_ta.e (000048 ) | § = exp (-.000057) = .99995 


[(.9997) (.9998))7" (.99995)° = .9850 
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1.4 RELIABILITY UNDER TEST CONDITIONS 
The ability of the system to pass the reliability test can be 
assessed by extrapolating the predicted reliability for a full 
system to provide an estimate of the MTBF for a reduced system 
under reliability test conditions. The reliability test will 
probably be performed on a system consisting of the following: | 
1 Position Computing Control 
1 Reference Position Set 
6 Positioning Sets 
3 DDU's 
The predicted reliabilities for test conditions are; 
PCC: Ri 000 = ,9985 


RPS: Rooo0 = ,9991 


PS R = ,9997 
(one): 3100 

R3200 = .9998 

R3300 = .9997 


The estimated reliability for the reduced system is: 
R, = (.9985) (.9991) '|(.9997) (.9998)|® (.9997)3 = .9928 
Assuming a 2-hour mission: 
.9928 = exp| -2 | 
2A= .0072, 
A = .0036 f/hour 


MTBF —e 280 hours 
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APPENDIX F-4 


RELIABILITY PREDICTION DATA SHEETS 


AN/PSQ-101 
AN/PSQ-101 
AN/PSQ-101 
AN/ASQ-148 
AN/ASQ-1 48 
AN/ASQ-148 


WORST CASE 

TYPICAL OPERATION 
RELIABILITY TEST CONDITIONS 
WORST CASE 

TYPICAL OPERATION 
RELIABILITY TEST CONDITIONS 
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APPENDIX F~4.1 

WORST CASE 
RELIABILITY PREDICTION DATA SHEETS 
POSITIONING SET AN/PSQ-101 
ENVIRONMENT: VEHICLE MOUNTED GROUND 
TEMPERATURE: +71°C 
STRESS LEVELS: DERATING LIMITS 
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Fermented 


APPENDIX F-4,2 


TYPICAL OPERATION 


RELIABILITY PREDICTION DATA SHEETS 
POSITIONING SET AN/PSQ-101 
ENVIRONMENT: VEHICLE MOUNTED GROUND 
TEMPERATURE: +50°C 

STRESS LEVELS: 30% (ASSUMED) 
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APPENDIX F-4.3 
RELIABILITY TEST CONDITION 


RELIABILITY PREDICTION DATA SHEETS 
POSITIONING SET AN/PSQ-101 
ENVIRONMENT: GROUND 

TEMPERATURE: +25°C 

STRESS LEVELS: 30% (ASSUMED) 
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pier: 


APPENDIX F-4.4 
WORST CASE 


RELIABILITY PREDICTION DATA SHEETS 

REFERENCE POSITION SET AN/ASQ-148 

ENVIRONMENT: AIRBORNE 

TEMPERATURE: +71°C 

STRESS LEVELS: DERATING LIMITS EXCEPT 
RUBIDIUM STANDARD STRESSES 
10% (ASSUMED) 
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APPENDIX F-4.5 
TYPICAL OPERATION 


RELIABILITY PREDICTION DATA SHEETS 
REFERENCE POSITION SET AN/ASQ-148 
ENVIRONMENT: AIRBORNE 

TEMPERATURE: +50°C 

STRESS LEVELS: 30% (ASSUMED) EXCEPT 


RUBIDIUM STANDARD STRESSES 
\ 


AT 10% (ASSUMED) 
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APPENDIX F-4.6 
RELIABILITY TEST CONDITION 


RELIABILITY PREDICTION DATA SHEETS 

REFERENCE POSITION SET AN/ASQ-148 

ENVIRONMENT: GROUND 

TEMPERATURE: +25°C 

STRESS LEVELS: 30% (ASSUMED) EXCEPT 
RUBIDUM STANDARD STRESSES 
10% (ASSUMED) 
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APPENDIX G 


MAINTAINABILITY PREDICTION DATA 


Vhs 

Data developed in performing the maintainability predictions 
are included in this appendix. This includes a complete set of 
maintenance task flow diagrams, task time summaries, and support- 
ing data. 
2. PREDICTION PROCEDURE 

The maintainability prediction was performed using a task 
flow analysis and time estimation technique involving the follow- 
ing steps: 

a. A maintenance task flow diagram is proposed using the 
symbology illustrated in Figure G-1l. The chart identifies 
each action performed and logical decision that must be 
made during each different maintenance action sequence. 
This diagram contains the following. 

(1) Identification of each action that consumes main- 
tenance time. 

(2. Identification of each logical decision and action 
taken following the decision. 

(3) failure rates of each module with respect to given 
failure symptoms and maintenance task sequences. 

b. Maintenance tasks are summarized on a worksheet and task 


—_— —s —— 0060 C06 ee 6 SU me ween 


time are summed to determine Mop values for each module 


and failure mode. These values are then weighted by the 
relative failure rate of the module for the respective 
failure mode, the weighted values summed, and the results 


divided by the total failure rate to obtain an M esti- 


CT 


mate. 
3. ASSUMPTIONS 5 
The maintainability prediction is based on the following 


assumptions: 


cin z No 
TASK TITLE aa 


LEGEND: 


Q. PREPARATION IR TEST TASK. BOK CONTAINS : 
TASK 1OENTIFICAT/ON NUMBER 
DESCRIPTIVE TASk TITLE 


b. DECISION. BOX INDICATES TYPE OF DECISION MADE AND NE&xr 
ACTION TAKEN FOLLOWING THE DECISION 


C. REPLACEMENT TASK, BOK CONTAINS: 
TASK IDENTIFICATION NUMBER 
IDENTIFICATION OF REPLACED ITEM 
REPLACED ITEM FAILURE RATE (2) 


d. RE-TEST ROVTINE. BOK INDICATES ACTION TAKEN FOLLOWING A 
REPLACEMENT. THE ABOVE EXAMPLE S4Y¥S: 


RE- PERFOR! TASK 3 
FAKE OK" DEC'SION INDICATED BY BLOCK h, 

GO TD TASK IF IF OR "DECISION IS REACHED 

GOT) SUBROUTINE ID IF "NOT OK” DECISION 1S PFALEED 


@ FLOW CHART ENTRY POINTS! ALABIC NUMERALS INDICATE A TASK 
IDENTIFICATION NUMBE?) ROMAN NUMERSIS INDICATE AN 
ENTRY PYINT IN ANOTHER SUGROYTINE 


4. THE DELAVS Nor ASSOC/ATED Will NUMBERED TASKS. B.G., TIME 
REQUIRED FOR TEST SEI 79 INDICATE AN EROIR ASTER 
PNE HAD BEN ENCOUNTERED 


Figure G-1. Maintenance Task Flow Diagram Symbology 


=e 


sepe oe 


Maintenance time starts after units are placed on the 
maintenance bench and all necessary equipment and spares 
are immediately available. 

All test equipment, including the TCS and auxiliary test 
equipment are turned on and ready for tests to begin. 
The time required to connect the unit to the test set is 
considered as part of the maintenance time. However, 
warm-up time, such as the time required for oscillator or 
frequency standard stabilization, is not considered. 

A maintenance action is completed when correction of the 
specific failure under consideration has been verified 
and the equipment returned to serviceable condition. 
Additional time that may be required to completely check 
out the remaining functions of the unit are not considered 
as part of a given corrective maintenance action. 

The following test equipment is assumed to be available; 
Test and Calibration Set. This provides facilities to 
test all digital circuitry of the DPU and the overall 
operation of the RTU. 

RTU Test Box. This provides test point "break-out" 
capability for using auxiliary test equipment. 
Oscilloscope. This is used for observing certain RTU 
signal status indications. 

Multimeter. This is used for measuring power supply and 
signal voltage levels. 


MAINTAINABILITY PREDICTION DATA 
Data prepared in performing the maintainability prediction are 
presented in the following sections of this appendix: 


Figure G-2. Positioning Unit Maintenance Flow 
Diagram 

Figure G-3. RTU Maintenance Flow Diagram (Receive 
Functions) 

Figure G-4. RTU Maintenance Flow Diagram (Transmit 
Functions) 

Figure G-5. Reference Position Unit Maintenance 


Flow Diagram 
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Figure G-6. Positioning Unit Supplementary 


Maintenance Flow Diagram (P7 Message 


Sequence) 
Table G-l. Corrective Maintenance Task Time 
Estimates 
Taole G-2. Positioning Set Task Time Summary 
Table G-3. Reference Position Set Task Time 
Summary (Defined Portions Only) 
Table G-4. RPS Fault Detection Time 
Table G-5. Elemental Task Times 


5S. REFERENCE POSITION SET MAINTAINABILITY ESTIMATE 

Table G-3 summarizes the maintainability prediction data for 
those portions of the RPS that have been sufficiently defined to 
date. However, the DPU has not been defined sufficiently to 
associate functions and failure modes with physical modu lization. 
Therefore, the predictions cannot be completed for these portions. 

An estimate of the maintenance time for the DPU is made by 
determining the approximate times expended before different 
categories of failures are detected by the TCS. In general, any 
DPU failure involves one of the following types of operation: 

a. Decoding and responding to commands. 

b. Generating and transmitting commands. 

c. Storing data in memory and reading data out of memory. 

Table G-4 indicates the approximate minimum, maximum and 
average times that are required to detect failures in each of 
these categories. This includes the time required for the TCS 
to lock up on an error that would be detected at the earliest 
possible and latest possible events of the test cycle. Assuming 
two modules would be replaced before the failure was corrected, 
the test would be performed three times for an "average'' DPU 
failure; twice for fault isolation and once for final check. 
Thus, the total troubleshooting time for a given categroy of 
failure is estimated to be the three times the average time 
indicated in Table G-4. 

It is assumed that there are 9 logic modules, each having a 


failure rate similar to those of the PS, or approximately 8 £/10°, 


aes 


e 


The memory failure rate has been predicted to be 98 £/10°, 
Assuming two modules are replaced on the average in locating a 
failure, and set-up time is 15.7 minutes (tasks 60, 61, and 64) 
the DPU maintenance time is calculated as: 


= 915.7 + 3(1.95)] (8) + [15.7 + 3(3.9)]98 
Mop 


24.8 minutes 


Combining this with the predicted values for the defined portions 


of the RPS (Table G-3) gives an estimate of Mop for the RPS as: 


31.2 x 639 + 24.8 x 170 — 99:18 winates 
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Figure G-2, Positioning Unit Maintenance 
Flow Diagram 
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G-13, G-14 


4NO2 2 ONINIL ‘8 
¥4009NZ 109UAS 'P 
willvWwYyos YLYO ‘2 
727$ 14AdNI AYOWIW * 
AYOUAW : 


g 
Vv 
JINAAGIS wi (BF 


TOYULNOD “SNVAL'P 
INO) ONY SNIWIL‘*D 
YIQOINA TOGWAS 'F 
YILLYWYOI WLlvd ‘oO 


1W99 SNvtl ‘9 
OD 3 Iwimi'y 


ee 


CMD Sw. TO 
ae 


52) press 
“SEND Msc." 


tags Jae, COMMAND 
aes \? CODER 


<2 ne 


emmieciiere nar oS 
J MES ACE ai 


OVIPYT “ONT 


: DDU O1sPLAy 


‘eerre 


Figure G-6. Positioning Unit Supplementary Maintenance 
Flow Diagram (P7 Message Sequence) 
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Table G-1. Corrective Maintenance Task Time Estimates 


Time 
Task Task Element (Min) 


Connect PU to Test Set, set TCU 
Switches 
Remove Unit from Case 


Latches (4) 

Captive Screws (4) 

Loose Screws (4) 

Multipin Connector (3) 

RF Connector (2) 

Pushbutton on Toggle Sw (12) 
Rotary Switch (3) 


Press "Send Message'' Press Switch +8 sec Hoa, 


Unplug PCB 
Exchange PCB 


Remove Protective Cover 
Connect to TCS 


Operate Switches 


07 Set Cond Switch to P4 Rotary Switch 0.2 | 
Press "Send Msg"' Press Switch +10 sec | 9.3 | 


rp bente ka Vigaet Gia apg p ee 


ery 


Table G-1. 


Check RTU Power 


Replace RTU Power Conv 


Set Xmit Cont to CONT 


Replace Rec-Trans Synth 


Corrective Maintenance Task Time Estimates (Cont) 


Cannon Connector 
Exchange Panel 
Set Switches (10) j 


Rotate Switch (6 positions) 
Meter measurement (6) 


Cannon Connector (2) 
Captive Screw (4) 
Exchange Module 


Flip Switch 
Wait 10 seconds 


Cannon Connector 
RF Connector (6) 
Captive Screw (2) 
Exchange Module 


Press Button 
Wait 10 seconds 


Rotate Switch (3 positions) 
Scope measurement (2) 
Mcter mcasurement 


Table G-1. Corrective Maintenance Task Time Estimates (Cont) 


ond Time 
Task Task Flemecnt (Min) 


F Replace RF/IF 


Replace Carricr IQ Loop 


Cannon Connector 
RF Connectors (5) 
Captive Screws (2) 
Exchange Module 


Cannon Connector 
RF Connectors (6) 
Captive Screws (2) 
Exchange Module 


Cannon Connector 
RF Connector (9) 
Captive Screws (4) 
Exchange Module 


Rotate Switch 
Meter measurement 


Cannon Conncctor 
Captive Screws (2) 
Exchange Module 


Replace CACU Sec Task 25 


Rotate Switch 
Meter measurement 


Sec Task 24 


Replace Code Detector 


Check Carrier Lock 


Replace Carrier 1Q Loop 


Table G-1. Corrective Maintenance Task Timc Estimatcs (Cont) i 


res eee 
No. Task Task Element (Min) 


ae | Replace CACU See Task 25 


Check Delayed Clock Lock 


Rotate Switch 
Meter measurement 


Cannon Connector 
RF Connectors (6) 
Captive Screws (2) 
Exchange Module 


a ee 
0 
1 


3 
Replace Code Detector See Task 27 
Replace CACU See Task 25 


42 Interrupt Transmit 


Replace RTU Cablic Harness Cannon Connectors (11) 


RF Conncctors (34) 


Pushbutton 
Wait 8 scconds 


Operate Self Test 


Table G-1. Corrective Maintenance Task Time Estimates (Cont) 


a Time 
Task Task Element re 


Check RTU Power See Task 16 
oo Replace RTU Power Connector See Task 17 


Replace PA & Modulator 


Cannon Connector 
RF Connector (3) 

Captive Screws (6) 
Exchange Module 


Total 


eC 
ee et 

Ca eae eae 
fe [ aimee | eer ia 
Ea Se 


Unplug-Plug PCB 
a PCB 


Replace Command Decoder 


= 


sa | Replace Message Output Cont Same as Task 53 a 
et Read DDU Display Digital Readout (10) a 


Replace DDU Exchange Unit 


Multipin Connector 


Connect RPS to Test Set, Set 
Switches 
Remove Units from Cases 


Captive Screws (8) 
Slide Unit out (2) (est) 


| ow om ow a aw le lel eles si ee 
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Table G~1. Correr, 2 Maintenance Task Time Estimates (Cont) 


Connect to TCS Multipin Connector (4) 
RF Connector (2) 
Push Button on Toggle (12) 


i Switch (3) 


Operate Switches 


i) =e 


61 Check DPU & Memory Voltages Meter measurements (5) Bal 
62 


Replace DPU Power Converter Unplug, Plug PCB 


Exchange PCB 


Replace Memory P.S. Captive Screws (4) 


Unplug-Plug Unit 
Exchange Unit 


Table G-2. Positioning Set Task Time Summary 


Command Decoder 
(Decode 81) 


Command Decoder 
(UUT Transmit on) 


Table G-2. Positioning Set Task Time Summary (Cont) 


Message Output Cont 
(UUT Transmit on) 


Message Output Cont 
(Encode S5) 


Word Control 
(Encode S5) 


3 
3 
1 
2 
1 
1 
1 
2 


cS 


Data Assemb 01-12 & 11 
(RC Error) 


= oO 
= bd 


alo 
Cy 


Memory 
(RC Error) 


Table G-2. Positioning Set Task Time Summary (Cont) 


Replaceable Item 


Control Panel 
Dele 


RTU Power Conv 
(Receive Problem) 


Rec-Trans Synth 
(Receive) 


RF/IF Amp 


35.1 


01- ge & 20 35,1 
8.9 
0.2 


44.2 


01- zi & 20 44,2 
2.8 
a 


Carrier IQ Loop 
(Code Search) 


CACU 
(Code Search) 


-_ 


CACU 
(Status Output) 


Table G-2. Positioning Set Task Time Summary (Cont) 
Replaceable Item 


Code Detector 
(Code Search) 


Carrier IQ Loop 
(Carrier Search) 


CACU 
(Carrier Search) 


01-31 25.7 
12.8 
0.3 
01-32 & 31 38.8 
5.9 
0.3 
16.5 
12.8 


r 
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co 
° 
oo 
iy 
o 
oo 
oo 
oo 


Code Detector 


li 
mm CO 


CACU 
(Status Output) 


01-29 


© wo 
oa 


Clock Phase Control 01-18 


(J) 
for) 


37 


CACU 01-38 
(Clock Adj) 39 
38 


Code Detector 01-39 & 38 
(Delay Circuit) 40 
38 
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Table G-2. Positioning Set Task Time Summary (Cont) 


CACU 
(Status Output) 


RTU Power Converter 
(Transmit Problem) 


1140 


PA and Modulator 


3790 


Rec/Trans Synth 
(Transmit) 


01-47 & 44 
48 


iw) 

Alo © © 
woe 
wooo 


> 

i 
ie) 
4 


CACU 
(Trans Code) 


01-48 & 44 
49 © 
44 


= 09 
wa 
Las 
on 


Totals 706 


i) 
(ex) 
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Table G-3. Reference Position Set Task Time Summary 


(Defined Portions Only) 


Replaceable Item Task # 


DPU Power Converter 


Memory Power Supply 


RTU Power Converter 
(Receive Problem) 


Rec-Trans Synth 
(Receive) 


RF/IF Amp 


Carrier IQ Loop 
(Code Search) 


oNrFNOF 
wore AD & & bd bY 


Table G-3. 


CACU (Status Output) 


Code Detector 
(Code Search) 


Carrier IQ Loop 
(Carrier Search) 


CACU (Status Output) 


Reference Position Set Task Time Summary 
(Defined Portions Only) (Cont) 


Table G-3, Reference Position Set Task Time Summary 
(Defined Portions Only) (Cont) 


Clock Phase Control 


CACU (Clock Adj) 


Code Detector 
(Delay Circuit) 


CACU (Status Output) 


RTU Power Converter 
(Transmit Problem) 


PA & Modulator 


Rec/Trans Synth 
(Transmit) 


WRT eo 


Table G-3. Reference Position Set Task Time Summary 
{Defined Portions Only) (Cont) 


CACU (Transmit 60-48 & 44 
Coder) 49 
44 


Table G-4, RPS Fault Detection Time 


Time After Starting Test 


0.1 3.0 


Decoding & Responding to Commands 


Generating & Transmitting Commands 


All Logic Modules 


Exercising Memory 


Table G-5. Elemental Task Times 


Fastener Operation 


NAVSHIPS 94324 
NAVSHIPS 94324 
NAVSHIPS 04324 
NAVSHIPS 94324 


Loose Screw with Washer 
DZUS Fasteners 

Captive Screw 

Latches 


Connector Operation 
pecs Li helio hedettabed t9 


Multipin Connector 

(Cannon Connector, etc.) 
(Plug-in) 

(2 Captive Screws) 


MIL-HDBK-472 
NAVSHIPS 94324 


eee eo ULC eae vs 


Table G-5, Elemental Task Times (Cont) 


RF Connector Disconnect 0.1 Motorola Timed Data 
Reconnect 0.4 


PCB 0.45 NAVSHIPS 94324 


Module Exchange 0. 20 


(Remove Module) 0.05 Estimate 
(Place On Chassis) 0.15 MIL-HDBK-472 


Task Element 1st Action Subsequent 
Test Operations 


Observe Go-No-Go 
(or digital readout) 


Estimated 


Monitor Analog Reading Estimated 
(Meter) 


Monitor Digital Reading Estimated 
(Scope) 


Pushbotton or 
Flip Switch 


Estimated 


Rotate Knob Estimated 
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